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_Jf y This  study  utiliiss^  Che  RELAP4/MOD5  Computer 

Code  to  analyze  the  effects  of  critical  flow 


modeling  on  the  thermal-hydraulic  transient 
response  of  a Ceneral  Electric  Soiling  water 
reactor  to  a major  primary  coolant  line  rupture. 
Included  TIT-  the  study  is  a presentation  of  the 
equations,  assumptions,  and  limitations  of  the 
critical  flow  models  available  for  use  in  REJAP4. 
Additionally,  an  evaluation  of  a temporary 
solution  to  a RELAP4  coding  error  associated  with 
stagnation  properties  calculation  is  presented. 

k • 

The  results  of  this  investigation  indicate 
that;  (l)  A solution  to  the  stagnation  properties 
calculational  error  that  exists  in  the  RE LAP 4 
code,  when  applied  to  the  Evaluation  “ode  1 , 
(RELAP4-EM)  provides  a conservative  evaluation, 
relative  to  the  Standard  Model  ( R ELA  Pu-SM ) , 
regardless  of  which  of  the  five  critical  flow 
models  available  in  the  code  is  selected;  (2)  Of 
the  five  critical  flow  monels  available  for  use  in 
RELAP4,  the  Moody  model  and  the  Hen r v -Fa us k e model 
are  nearly  equivalent  in  their  relative  degrees  of 
conservatism.  Following  these  models,  in  order  of 
relative  degree  of  conservatism,  are  the  “od  If  led 
Momentum/Homogeneous  Equilibrium  model,  the 
Homogeneous  Equilibrium  model,  and  the  Sonic 
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In  order  to  best 


alleviate  the 


stagnation  properties  coding  error,  the  flow  areas 
for  the  volumes  Immediately  upstream  from  the 
recirculation  line  rupture  should  he  increased 
such  that  the  area  ratio  (break  flow  area  / 
upstream  volume  flow  area)  is  0.7. 
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ABSTRACT 

This  study  utilizes  the  RELAP4/MOD5  Computer 
Code  to  anslyze  the  effects  of  critical  flow 
modeling  on  the  thermal-hydraulic  transient 
response  of  a General  Electric  boiling  water 
reactor  to  a major  primary  coolant  line  rupture. 
Included  in  the  study  is  a presentation  of  the 
equations,  assumptions,  and  limitations  of  the 
critical  flow  models  available  for  use  in  RELAP4. 
Additionally,  an  evaluation  of  a temporary 
solution  to  a RELAP4  coding  error  associated  with 
stagnation  properties  calculation  is  presented. 

The  results  of  this  investigation  indicate 
that;  (1)  A solution  to  the  stagnation  properties 
calculations  1 error  that  exists  in  the  R FLAP  4 
code,  when  applied  to  the  Evaluation  Model, 
(RELAP4-EM)  provides  a conservative  evaluation, 
relative  to  the  Standard  Model  (RELAP4-SM), 
regardless  of  which  of  the  five  critical  flow 
models  available  in  the  code  is  selected;  (2)  Of 
the  five  critical  flow  models  available  for  use  in 
RELAP4,  the  Moody  model  and  the  Henry-Fauske  model 
are  nearly  equivalent  in  their  relative  degrees  of 
conservatism.  Following  these  models,  in  order  of 
relative  degree  of  conservatism,  are  the  Modified 
Momentum/Ho mo geneous  Equilibrium  model,  the 
Homogeneous  Equilibrium  model,  and  the  Sonic 
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model;  (3)  In  order  to  best  alleviate  the 
stagnation  properties  coding  error,  the  flow  areas 
for  the  volumes  Immediately  upstream  from  the 
recirculation  line  rupture  should  he  Increased 
such  that  the  area  ratio  (break  flow  area  / 
upstream  volume  flow  area)  Is  0.7, 
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1.  INTRODUCTION 


1 . 1 GENERAL 


The  United  States,  in  its  fourth  decade  of  nuclear 
reactor  research  and  operation,  now  has  seventy-two  civilian 
power  reactors  in  operation  and  has  generated  over  a 
trillion  megawatt-hours  of  electricity.  Since  the  Inception 
of  the  nuclear  power  industry,  safety  has  been  the  major 
consideration  in  the  design,  construction,  and  operation  of 
nuclear  power  plants. 

Modern  experimental  programs  in  water-reactor  safety 
are  directed  toward  improving  the  ability  to  understand 
specific  phenomena  and  the  results  of  tests,  with  the  goal 
of  applying  this  understanding  to  the  interpretation  of 
"similar"  phenomena  that  might  take  place  if  a reactor 
accident  were  to  occur.  Complex  digital  computer  codes  have 
been  developed  to  predict  the  t he rma 1 -h yd r a u 1 ic  transient 
response  of  nuclear  reactor  systems  to  hypothesized 
accidents.  Verification  of  specific  portions  of  these 
computer  codes  is  the  goal  of  an  ongoing  experimental  water- 
reactor  safety  program. 


1.2  LOSS-OF-COOLANT  ACCIDENT  (LOCA) 


A LOCA  Is  generally  defined  as  any  accident  which 
causes  the  coolant  to  be  lost  from  the  cooling  system 
serving  the  reactor  core.  The  design-basis  LOCA  Is  that 
primary  system  pipe  rupture  which  results  In  maximizing  the 
fuel-rod  cladding  temperatures  attained  during  the  accident. 

For  boiling  water  reactors  (BWRs)  the  design-basis  LOCA 
is  a double-ended  guillotine  break  of  one  of  the  two  coolant 
recirculation  lines.  The  coolant  recirculation  lines  are 
the  primary  means  of  maintaining  coolant  flow  through  the 
reactor  core,  and  are  the  largest  diameter  piping  in  the 
primary  system.  Emergency  core  cooling  systems  tECCS)  are 
designed  to  provide  cooling  to  the  core  in  order  to  prevent 
fuel  element  melt  down  In  the  event  of  a design-basis  LOCA. 
Future  references  to  LOCA  will  Imply  design-basis  LOCA 
unless  otherwise  indicated. 

Within  the  present  regulatory  framework,  the  thermal- 
hydraulics  of  a LOCA  can  be  described  as  having  three 
distinct  phases  known  as  blowdown,  refill,  and  reflood  (2). 
(1)  Reactor  Blowdown:  During  the  reactor  depressurization  or 
blowdown  phase,  reduction  of  core  flow  and  increase  in  the 
local  quality  result  in  departure  from  nucleate  boiling 
(DNB)  and  dry-out  of  the  core.  Following  the  establishment 
of  the  conditions  causing  DNB,  the  heat  transferred  from  the 
fuel  rods  to  the  coolant  decreases  markedly,  and  the 


3 


temperature  of  the  fuel  rods,  and  the  fuel-rod  cladding, 
rises  rapidly.  Depressurization  continues  until  primary 
system  and  containment  pressures  are  equalized  and  the  break 
flow  Is  essentially  zero. 

(2)  Refill  and  Reflood:  At  the  end  of  the  blowdown  phase, 
the  liquid  coolant  Inventory  In  the  system  Is  small,  and  the 
core  may  be  completely  immersed  In  steam.  The  system 
pressure  Is  low,  on  the  order  of  50  to  100  PSIA.  As  the 
ECCS  starts  injecting  emergency  coolant,  refilling  of  the 
lower  plenum  begins.  For  conservative  safety  evaluations, 
heat  transfer  to  the  coolant  Is  often  totally  Ignored  during 
this  refill  phase.  When  the  liquid  level  reaches  the  bottom 
of  the  core,  the  reflooding  phase  begins.  As  the  liquid 
level  rises  In  the  core  a tremendous  amount  of  heat  Is 
removed  from  the  fuel  rods  due  to  vaporization  of  the 
Incoming  liquid.  Peak  cladding  temperatures  will  he 
attained  during  this  phase  of  the  accident. 

For  an  excellent,  detailed  physical  description  of  a 
BWR  LOCA  the  Interested  reader  Is  referred  to  references  (3) 
and  ( 4 ) . 

1 . 3 CODE  DESCRIPTION 

RELAP4  is  a computer  program,  written  in  Fortran  IV, 
that  was  developed  primarily  to  describe  the  thermal- 
hydraulic  transient  behavior  of  water-cooled  nuclear 
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reactors  subjected  to  postulated  accidents  such  as  those 
resulting  from  loss  of  coolant,  pump  failure  or  nuclear 
power  excursions.  Fundamental  assumptions  Inherent  In  the 
thermal-hydraulic  equations  used  In  the  code  are  that  a two- 
phase  fluid  Is  homogeneous  and  that  the  phases  are  In 
thermal  equilibrium.  Models  are  available  to  modify  the 
homogeneous  assumption. 

The  program  Is  currently  available  in  Its  fourth  major 
revision,  hence  RELAP4.  The  edition  of  the  program 
Installed  at  The  Pennsylvania  State  University  Is 
RELAP4/MOD5-UPDATE2  (henceforth  referred  to  as  RELAP4). 

There  are  four  major  options  available  In  RELAP4.  The 
Evaluation  Model  option,  hereafter  referred  to  as  RELAP4-FM, 
Is  Intended  to  provide  conservative  results  which  meet  most 
of  the  requirements  for  licencing  under  Appendix  K to  10  CPR 
Part  50.  Another  option  available,  called  the  Standard 
Model  (RELAP4-SM),  Is  used  to  provide  the  best  estimate  of 
the  actual  expected  response  during  a LOCA.  While  RELAP4 
was  primarily  designed  to  be  a blowdown  code,  it  Is  capable 
of  representing  the  LOCA  transient  through  reflood  In 
pressurized  water  reactors  (R  ELAP4-FLOOD  option).  Also 
available  Is  the  containment  option  (R ELA P4-C0NTA INMENT ) , 
which  can  be  used  to  evaluate  the  effectiveness  of  Ice- 
condenser  containment  systems.  The  two  main  options  that 
will  be  used  in  this  study  are  RELAP4-EM  and  RELAP4-SM. 
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RELAP4  requires  numerical  Input  data  which  completely 
describes  the  Initial  conditions  and  the  geometry  of  the 
system  being  analyzed*  Included  In  the  Input  data  are 
physical  characteristics  such  as  fluid  volume,  geometry, 
pump  characteristics,  power  generation,  and  material 
composition*  Transients  can  be  Initiated  by  the  control 
action  Inputs  to  the  program*  These  controls  can  describe 
breaks  in  fluid  piping,  core  power  level  variations 
Including  a SCRAM,  valve  actions,  etc.  For  each  time 
advancement  the  program  computes  fluid  conditions  such  as 
flow,  pressure,  quality,  and  mass  Inventory*  Also  computed 
are  thermal  conditions  within  the  solid  materials  such  as 
temperature  profiles  and  power  generated,  and  the  fluid- 
solid  Interface  conditions  such  as  heat  flux  and  surface 
temperature  * 

A complete  theoretical  development  of  the  code,  as  well 
as  implementation  procedures,  Is  contained  in  the 
R EL A P 4 /MOD  5 Users  Manual  (5).  An  excellent  description  with 
specific  reference  to  BWR  LOCA  analysis  Is  also  available  In 
a masters  thesis  by  Laird  (3). 

The  reactor  system  used  for  this  study  Is  the  Hope 
Creek  Nuclear  Power  Generating  Station*  This  General 
Electric  BWR  was  modeled  for  RELAP4  by  Idaho  National 
Engineering  Laboratory  (INEL)  personnel,  and  was  Included  as 
a sample  problem  with  the  RELAP4  transmittal  tape. 
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The  thermal-hydraulic  model  of  Hope  Creek,  depicted  In 
Figure  1,  consists  of  twenty-three  control  volumes,  thirty- 
three  Junctions,  and  twenty-five  heat  slabs.  Table  1,  pages 
8 and  9,  presents  a description  of  the  twenty-three  control 
volumes  used  In  this  model.  Junctions  simply  connect 
control  volumes,  and  heat  slabs  represent  sections  of  the 
physical  system  which  can  absorb,  radiate,  store,  or 
generate  heat.  RELAP4-EM  and  RELAP4-SM  data  Input  decks 
used  In  this  study  are  discussed  in  Appendix  A. 


1.4  CRITICAL  FLOW 

"When  a fluid  expands  from  a compressed  state  to 
arbitrary  ambient  (or  "receiver")  conditions  in  passage 
through  an  outlet,  the  flow  rate  Is  always  less  than  a 
certain  maximum  or  "critical"  value.  After  this  critical 
condition  has  been  reached,  further  reductions  In  receiver 
pressure  leave  the  flow  rate  unaltered,  serving  only  to  form 
a steep  pressure  gradient  at  some  location  in  the  outlet 
passage  (6)."  This  condition  Is  commonly  called  critical, 
or  choked,  flow. 

At  the  Instant  of  recirculation  line  rupture  the  flow 
rate  out  of  the  break  Increases  until  critical  flow 
conditions  are  reached.  There  are  five  critical  flow  models 


available  In  RELAP4  to  predict  this  critical  flow  rate. 
They  are,  (1)  The  Sonic  model,  (2)  The  Moody  model,  (3)  The 


FIGURE  1 

RELAP4  MODEL  OF  G.  E.  BWR 
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Table  1.  Volume  Descr iptions  in  RE LAP 4 Model 
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Volume  No.  Description 


1 Upper  Plenum  and  Steam  Separators 

2 Upper  Downcomer  (This  volume  extends  from  the  feedwater 
entry  point  to  the  bottom  of  steam  dryers.) 

3 Steam  Dome  (This  volume  contains  the  high  quality  steam 
located  above  the  dryers . ) 

4 Lower  Downcomer  (This  volume  contains  the  subcooled  water 
which  enters  the  downcomer  from  the  feedwater  entry 
point.  The  lower  downcomer  describes  the  region  of  the 
downcomer  between  the  feedwater  entry  point  and  the 
entry  point  of  the  feedwater  into  the  lower  plenum.) 

5 Broken  Loop  Recirculation  Suction  Line  (This  volume 
models  the  26"  pipe  which  pulls  the  saturated  water 
from  the  lower  downcomer  into  the  jet  pump  of  the 
broken  loop.  Although  the  break  occurs  along  this  suc- 
tion line,  the  actual  break  volume  is  modeled  by 
volume  21.) 

6 Broken  Loop  Recirculation  Pump 

7 Broken  Loop  Recirculation  Discharge  Line  (This  volume 

models  the  24"  line  which  returns  the  water  from  the 
recirculation  pump  to  the  jet  pump  intake.) 

8 Intact  Loop  Recirculation  Suction  Line 

9 Intact  Loop  Recirculation  Pump 

10  Intact  Loop  Recirculation  Discharge  Line 

11  Lower  Plenum 

12-16  Heated  Core  Sections  (With  the  addition  of  volume  23, 

these  volumes  model  the  core.  To  achieve  relatively 
equal  thermodynamic  characteristics,  the  heated  core 
sections  vary  in  height  from  2.8S  feet  to  1.00  feet.) 

17  Guide  Tubes  (The  control  rod  guide  tubes  in  the  lower 

plenum  and  the  unheated  segment  of  the  core  are  modeled 
as  one  control  volume.) 
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Table  1 . (cunt. ) 


Vo  1 ume  No . 

Description 

Id 

Core  Bypass  (This  volume  represents  the  region  between 
the  fuel  assemblies  in  the  core.) 

19 

Intact  Loop  Jet  Pump  (This  volume  models  the  intact  loop 
]et  pumps  from  the  intake  to  the  exhaust  into  the  lower 
downcomer.  Each  bank  of  ten  jet  pumps  is  modeled  as 
one  volume . ) 

20 

Broken  Loop  Jet  Pump  (This  single  volume  models  the 
bauik  of  ten  jet  pumps  in  the  broken  loop.  ) 

21 

Break  Volume  (This  volume  is  specified  between  the 
broken  loop  recirculation  discharge  line  and  the  lower 
downcomer. ) 

22 

Containment 

23 

Unheated  Segment  of  the  Core 
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Henry-Fauske  model,  (4)  The  Homogeneous  Equilibrium  model, 
and,  (S)  The  Modified  Momentum/Homogeneous  Equilibrium 
model.  These  models  are  discussed  In  detail  in  Chapter  Two 
of  this  study. 

The  simple  momentum  equation  for  inertial  flow  rate  is 
quite  accurate  at  low  flow  rates  but  becomes  greatly 
exaggerated  as  the  driving  pressure  differential  Increases 
and  critical  conditions  are  approached.  RELAP4  computes  the 
flow  rate  using  the  simple  momentum  equation  and  compares 
this  flow  rate  to  that  predicted  by  the  selected  critical 
flow  model.  The  minimum  is  then  selected  as  depicted  by  the 
solid  line  in  Figure  2 (7). 


FIGURE  2 

JUNCTION  MASS  FLOW  RATE  SELECTION 
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1.5  STAGNATION  PROPERTIES 

All  of  the  critical  flow  models  available  In  RELAP4, 
with  the  exception  of  the  Sonic  model,  require  that 
stagnation  properties  be  calculated  for  the  volume  upstream 
of  the  Junction  undergoing  critical  flow.  These 
calculations  are  performed  by  a subroutine  called  STAGP. 

The  stagnation  pressure  and  enthalpy  are  calculated 
from  the  static  pressure  and  enthalpy  using  a calculated 
fluid  velocity.  An  isentropic  path  from  the  static  to  the 
stagnation  enthalpy-pressure  point  Is  followed. 

The  stagnation  enthalpy  is  calculated  from  the  kinetic 
energy  relationship: 

hQ  - hL  + V2/2gcJ  (1) 

where  V Is  calculated  from  the  following  relationship: 

V • Wv/A  (2  ) 

The  quantities  h^,  W,  and  v are  evaluated  at  the  Junction,  A 
Is  the  upstream  volume  flow  area. 

In  equation  (1),  V is  not  allowed  to  he  greater  than 
the  Isentropic  sonic  velocity  of  the  upstream  volume. 

Applying  the  first  law  of  thermodynamics  In 
differential  form  for  a fluid  undergoing  a flow  change 
provides,  (Gibbs  equation); 


Tds  « dh  - vdP  . 


For  an  isentroplc  process,  ds  - 0.  Then 


d P - dh / v 


■T> 
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(3) 
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Integrating  this  for  the  change  In  enthalpy  and 
applying  static  and  stagnation  limits  provides: 

p.  ■ pi  + 3,°  <lh/v  (5> 

where  v,  the  specific  volume,  Is  a function  of  P and  h,  and 
Pq  and  Pj  are  the  stagnation  and  static  pressures, 
respectively. 

Equation  (5)  is  advanced  from  (h  ,P  ) to  (h  ,P  ) by  a 

11  o o 

fourth-order  Runge-Kutta  subroutine. 

It  had  been  reported  that  there  existed  a logic  problem 
in  the  stagnation  properties  subroutine,  STACP,  requiring 
the  elimination  of  this  routine  from  the  RELAP4  logic  (8, 
9).  It  has  since  been  discovered,  and  disclosed  at  the  INEL 
Code  Users  Workshop  conducted  at  INEL,  Idaho  Falls,  Idaho 
8/28/78  thru  9/1/78,  that  the  stagnation  properties 
calculation  problem  was  not  due  to  a logic  error  in  the 
STAGP  subroutine,  hut  rather  was  due  to  the  fact  that  the 
code  was  predicting  velocities  greater  than  the  isentroplc 
sonic  velocity  for  the  volume  immediately  upstream  of  the 
break  Junction. 


As  pointed  out  earlier,  the  code  does  not  allow  the 
calculated  velocity  In  equation  (1)  to  exceed  the  isentropic 
sonic  velocity  of  the  upstream  volume.  If  the  Junction 
conditions  are  such  that  equation  (2)  predicts  a velocity 
greater  than  sonic,  erroneous  values  for  the  stagnation 
pressure  and  enthalpy  will  result. 

The  "fix"  suggested  at  the  workshop  was  to  artificially 
Increase  the  upstream  volume  flow  area  (the  denominator  of 
equation  (2))  in  order  to  prevent  supersonic  velocity 
predictions  from  occurring.  This  flow  area  Increase  should 
be  carried  out  for  all  volumes  where  critical  flow  Is  likely 
to  occur. 

1.6  PURPOSE 


The  purpose  of  this  study  Is  threefold; 

(1)  Review  the  derivation  of  the  equations  for  each  of 
the  critical  flow  model  options  available  in  RELAP4. 

(2)  Determine  by  what  amount  the  upstream  volume  flow 
area  must  be  increased  in  order  to  alleviate  the  problem  in 
calculating  stagnation  properties  that  was  reported  in  the 
previous  section. 

(3)  Present  a sensitivity  study  on  the  critical  flow 
models  available  in  RELAP4  and  rate  the  models  as  to  their 
relative  degrees  of  conservatism. 


2.  CRITICAL  FLOW  MODELS 


2.1  INTRODUCTION 

As  previously  reported,  there  are  five  critical  flow 
models  available  in  RELAP4.  They  are: 

(1)  The  Sonic  model  (10), 

( 2 ) The  Moody  model  (11), 

(3)  The  Henry-Fauske  model  (12), 

(4)  The  Homogeneous  Equilibrium  model  (13),  and 

(5)  The  Modified  Momentum/Homogeneous  Equilibrium  model 
(14)  . 

For  each  junction  in  the  system  the  user  must  specify 
whether  the  code  is  to  use  one  of  these  models  or  the  simple 
momentum  (inertial  flow)  equations  to  estimate  the  flow  rate 
through  that  Junction.  In  addition  the  user  may  select  f ron 
eighteen  possible  combinations  of  these  models.  For  example 
the  user  may  elect  to  use  the  Sonic  model  when  the  fluid  is 
subcooled  and  the  Moody  model  when  the  fluid  is  saturated. 
These  combinations , as  well  as  the  selection  methods,  are 
discussed  in  the  RELAP4/MOD5  Users  Manual  (15). 


The 

purpose  of  this  chapter  is 

to  present 

each 

of  the 

critical 

flow 

models 

listed  above  , 

provid  ing 

the 

basic 

assumptions  and 

limitations  of  each. 

In  addition, 

at  the 

end  of 

this 

chapter 

a brief  description  of 

the 

logic 

involved  in  junction  flow  rate  estimation  is  given. 


2.2  THE  SONIC  MODEL  (10) 
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As  the  name  Implies,  this  model  assumes  that  the 
Junction  flow  rate  Is  that  which  will  result  In  sonic 
velocity  at  the  Junction. 

The  sonic  velocity  of  the  fluid  is  calculated  from, 

a2  - [Xv  +(l-X)v,]2g  T(P')2/{X[C_  J+TP'v  (P'K  -23  )] 
g f °c  P g a g 

8 

(6) 

+ ( 1 -X ) [ C p J + TP'vf (P'Kf  - 28  f ) J > 

where  P'  signifies  the  derivative  of  the  partial  pressure 
with  respect  to  temperature. 

It  Is  important  to  note  that  the  independent  variables 
in  equation  (6)  are  based  upon  the  upstream  volume 
conditions.  That  19,  the  steam  tables  are  entered  with 
upstream  volume  specific  internal  energy  and  specific  volume 
as  the  independent  variables.  All  other  upstream  volume 
parameters,  including  pressure  and  sonic  velocity,  are 
dependent  on  these  two  quantities.  The  result  is  that  the 
upstream  volume  critical  velocity  is  calculated  and  then 
impressed  upon  the  local  downstream  Junction. 

The  critical  mass  flow  rate  through  the  junction  is 
calculated  from  the  one-dimensional  continuity  equation, 

W • 0 A a ( 7 ) 


The  density  used  In  equation  (7)  Is  calculated  from  the 
upstream  volume  density  as  follows.  The  upstream  volume 
density  Is  first  adjusted  for  the  frictional  loss  and 
kinetic  energy  changes  to  the  Junction.  The  fluid  is  then 
expanded  1 se n t r op  1 c a 1 ly  from  the  Mach  number  Just  upstream 
from  the  Junction  to  that  In  the  Junction,  which  Is  Mach 
one,  when  sonically  choked.  It  should  be  noted  that  A^  Is 
the  flow  area  at  the  junction. 

The  assumptions  Implicit  In  the  use  of  this  model  are: 

(1)  The  critical  flow  rate  Is  assumed  to  be  that  which 
will  result  In  sonic  velocity  at  the  Junction, 

(2)  The  critical  flow  rate  Is  given  by  the  one- 
dimensional  continuity  equation, 

(3)  The  Junction  critical  velocity  is  based  upon 
upstream  volume  fluid  conditions,  and 

(4)  In  calculating  the  density  used  in  eouatlon  (7) 
isentroplc  expansion  Is  assumed. 

This  model  is  also  capable  of  predicting  critical  flow 
rates  when  the  upstream  volume  contains  air  or  a mixture  of 
air  and  saturated  or  superheated  steam.  This  capability  Is 
relevant  only  to  the  containment  option  ( R ELA  P4 -CONTA IVMFNT ) 


and  will  not  he  further  discussed 


1 7 


2 . 3 THE  MOODY  MODEL  (U  ) 

The  Moody  model  la  applicable  only  In  the  saturated 
region,  and  therefore  la  used  primarily  In  low  pressure 
systems.  It  is  the  only  critical  flow  model  available  In 
RELAP4  which  allows  for  slip  between  phases.  In  addition  it 
does  not  rely  on  empirical  data  for  Its  derivation. 

The  derivation  of  the  equation  for  critical  mass  flux 
starts  with  the  continuity  and  energy  equations  for  two- 
phase,  annular  flow; 

G - W/A  ■ aV  /Xv  ■ (1  - a ) V / ( 1 - X)v,  (P) 

a g f f 

h - X(h  + V2  /2g  J)  + (1  - X ) (h  + V2 /2  g J)  . (4) 

o g R C f f C 

X and  a are  defined  by, 


X - V /V  ; 1 - X - W./W  (10) 

P f 


■J  * A /A  ; 1 - a - a /A  (11) 

8 ^ 


The  ratio  of  average  vapor  velocity  to  average  liquid 
velocity  Is  the  slip  ratio,  S ; 


S - V /Vt 

K f 


(12) 


j 


2 
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Rewriting  equation  (8)  and  substituting  equation  (12) 
yields  , 

V /V , - (1  - a ) Xv  /c(l  - X ) v , - S . (13) 

g f g r 

Solving  for  the  void  fraction  in  terms  of  the  slip 
ratio  yields, 

a - 1/(1  + [ S v . ( I - X)/Xv  ) ) . (14) 

f g 

Assuming  that  the  expansion  is  isentropic  (s  - s^),  an^ 
from  the  definition  of  quality, 

X - (s  - s,  ) /s,  (15) 

o f f g 


o r 


1 - X ■ ( s - s ) / s - • (18) 

g 0 E g 

An  expression  for  the  mass  flux  Cl  is  obtained  by 
combining  equations  (8)  through  (16)  as  follows; 

Rewriting  equation  (9)  yields, 

2 g J (h  - h - X(h  - h ) 1 - V2[X  + V2(l-X)/V2) 

c o f g f g f g 


(17) 
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Substituting  equation  (12)  Into  equation  (17)  and 


noting  that  h ■ h - h,  yields, 
‘ * R f 


2 g J (h  - h.  - Xh  ] - V2[X  (1  - X ) / S 2 
cot  tg  g 


(18) 


Rewriting  equation  (14)  aa. 


m - Xv  / (Xv  + Sv  ( 1 - X)] 
R g f 


(19) 


solving  equation  (8)  for  V , and  substituting  equation  (l3) 

R 


yields  , 


V - C [Xv  + Svc(  1 - X) ] 
g g ‘ 


(20) 


Substituting  equation  (20)  Into  equation  (18)  and 
solving  for  G yields. 


G “ 2*cJ(ho  - hf  ’ '^fg^ 


(2  1 ) 


{ [Xv  + S v . ( 1 - X ) ) 2 [ X + (1  - X ) / S 2 J ) 
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Substituting  equations  (15)  and  (18)  into  equation  (21) 
y 1 e l i s , 


"1 


?n 


C - 2*cJIho  - hf  - hfg(So  - “f^fg1' 


( [Sv  (a  - a )/s  + v (s  - a )/s  ] 

f g o fg  go  f fg 


(22) 


((a  - a )/s  + (a  - a ) /S  s )) 

o ffg  g o fg 


Note  from  equation  (22)  that  If  h and  s are  known, 

oo 

then  0 i a a function  of  S and  P (ie.  knowing  the  pressure 

and  the  fact  that  saturation  exists  fixes  all  of  the 

quantities  in  equation  (22)  except  h , s , and  S). 

o o 

Therefore,  If  G has  a maximum  It  must  satisfy  the 
necessary  conditions, 


3C  n 
93  P " ° 


(23) 


and. 


9C 

7? 


(24  ) 


assuming  that  S and  P are  Independent* 

Applying  equation  (23)  to  equation  (22)  and  solving  for 
S yields , 

S • (v  /v  )1/3  . (25) 

f?  f 

Thus  for  maximum  flow  rate  S depends  only  on  P. 


When  equations  (25)  and  (22)  are  combined,  a maximum  C 
must  satisfy  equation  (24).  Equations  (22),  (25),  and  the 
state  equation. 


h 

o 


+ h 


f8, 


(s 


sf  ) /s 
0 


f*o 


(26) 


were  programmed  for  saturated  steam/water  properties  (16) 

and  solved  by  digital  computer  for  values  of  pressure  which 

satisfied  equation  (24)  for  Input  values  of  P and  h . 

o o 

Calculations  showed  that  G has  a single  maximum  value  for 
known  Pq  and  hQ.  The  maximum  mass  flux  was  thereby 
determined  . 

The  results  of  the  above  computations  are  available  In 
tabular  form  for  use  In  RELAP4.  The  tables  are  entered  with 
upstream  volume  stagnation  pressure  and  junction  stagnation 
enthalpy  and  a single  mass  flux  is  returned. 

The  Moody  model  Is  based  upon  an  annular,  two-phase, 
o ne -d  linen  s 1 o n&  1 model  with  uniform  axial  velocity  of  each 
phase  and  thermodynamic  equilibrium  between  phases. 
Additional  assumptions  are: 

(1)  Both  phases  experience  the  same  local  static 
pressure, 

(2)  The  flow  is  isentroplc  from  entrance  to  exit,  hence 
stagnation  enthalpy  Is  constant, 

(3)  The  liquid  phase  Is  Incompressible,  and 
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. 


"i 
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(4)  The  slip  ratio  and  pressure  are  independent 
variables. 

The  Moody  model  is  strictly  applicable  only  in  the 
saturated  regime  with  stagnation  pressure  limits  of  1 PSIA 
to  3000  PSIA.  If  the  Moody  model  is  used  in  the  suhcooled 
regime  the  code  will  enter  the  critical  flow  tables  at  zero 
quality. 


2.4  THE  HENRY-  FAUSKE  MODEL  (12) 

The  Henry-Fauske  model  was  developed  primarily  to 
account  for  the  nonequilibriura  nature  of  fluid  flow  under 
critical  conditions  and  is  applicable  for  fluids  in  the 
subcooled  as  well  as  the  saturated  and  superheated  states. 

The  derivation  starts  with  the  steady-state,  one- 
dimensional  continuity  and  momentum  equations  for  two-phase 
flow; 

Liquid  continuity. 


WfVf 


- Vf 


(2  7) 


Vapor  continuity. 


W v ■ A V (2f>) 

g R g g 
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Momen  tun , 

-AdP  - d(V  V + WfV()/n  +d  F . ( 2 P ) 

R g ft  c w 

For  high  velocity  flows  the  wall  shear  forces  are 

negligible  compared  to  the  momentum  and  pressure  gradient 

terms.  Therefore,  with  the  substitution  of  the  liquid  and 

vapor  flow  rate  definitions,  (W  - X W and  U - (1  - X)W) 

R f 

equation  (29)  becomes, 

-dP  - Cd[XV  + (l  - X)Vf]/g  . (30) 

g f e 

It  is  assumed,  for  fixed  upstream  stagnation 

conditions,  that  X,  v , v , V , and  V are  either  constant 

8 f R f 

or  composite  functions  of  P and  z.  Hence,  equation  (30) 
becomes  , 

Cj 1 - + (1  - x>vf]/gc>j  (31) 

where  the  subscript  J Indicates  that  the  enclosed  quantities 
are  evaluated  at  the  junction.  At  critical  flow,  the  mass 
flow  rate  exhibits  a maximum  with  respect  to  junction 
pressure  ; 


dCc/dP|j 


0 


(32) 


Equation  (32)  can  be  applied  to  equation  (31)  to  give 
an  expression  for  the  critical  flow  race; 

-G'2  - (jtsUXS  + (l-X)l  [ (l-X)Sv{  + Xvj?])/gcS>J  (33) 

where,  as  before,  S is  defined  by  S ■ V /V, . By  taking  the 

g f 

Implied  derivative,  equation  (33)  can  he  expanded  to 

-C'2  - <X[l+X(S-l)]dv  /dP  + S [1  + X (S-2 )-X2  (S-l ) ] dve /dP 
c g f 

+ v [l  + 2 X ( S - 1 ) + 2 S v (X-l ) + S ( l -2 X ) ] dX  /d  P (34) 
g ‘ 

+ X ( 1 -X ) [ v S - V /SJdS/dP)  /g  s . 

f 8 J c 

Equation  (34)  can  be  simplified  by  using  the  following 
assumptions  and  relationships: 

(1)  The  vapor  and  liquid  velocities  are  equal  (no 
slip); 

S - 1 (3M 


(2)  The  liquid  phase  Is  Incompressible; 
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(3)  The  vapor  behavior  can  be  described  by  a polytropic 
process  such  that, 

dv/dP*v/nP  (37) 

g g 

where  n Is  the  thermal  equilibrium  polytropic  exponent 

derived  by  Tangren  et  al.  (17)  and  given  by, 

n - [(l  - X)C„  /C„  + X]/[(l  - X)C„  /C..  + X/Y  ] (37a) 

Vj  v V g 

(4)  Negligible  interphase  mass  transfer; 

X - X ()«) 

J ° 

(5)  The  empirical  result; 

dS/dP|  - 0 . (3<>) 

Equations  (35)  through  (39)  serve  to  reduce  equation 

(34)  to 

G.  - [X  v /nP  - (v  - v,  )dX/dP]~l/2g  1/2  (40) 

Cog  g fQ  J C 

Formulating  the  actual  junction  fluid  quality  in  terms 
of  the  equilibrium  quality  at  that  location  and  making  use 
of  empirical  data  as  well  as  some  mathematical  manipulation, 
equation  (40)  becomes 


26 


1/2 


G - gWi  { X v /nP  + 
c c o g 


(v  - v.  ) [{  (1  - X )N/s.  >ds  /dP 

8 f„  ° f8p  fP 


(4  1 ) 


{X  C _ / Ps.  )(l/n  - 1/Y)]  > 71/2 

0 P„  J 

g o 


The  nonequilibrium  parameter  N is  of  primary  import  in 
the  low  quality  regime.  Rased  on  experimental  data  (IP), 
this  parameter  can  be  related  to  the  equilibrium  quality 


(X  ) as  follows , 

J 


N(XE  ) - 
'J 


X /0  . 1 4 , X <0.14 

J EJ 


(42) 


, XE  >0.14 


J 


The  relation  for  critical  mass  flux,  equation  (41),  is 
coupled  with  the  two-phase  momentun  equation; 


(1-X  ) v _ P ( 1 — R ) + Y X P (v  - Rv  )/(  Y-  1) 

o f o o o g g . 

o o J 


(43) 


c!|Vg  * vf0(1  ' V,2/2*e 


Solving  equations  (41)  through  (43)  simultaneously 
generates  the  Henry-Fauske  tables  used  in  FELAP4.  The 
tables  are  accessed  with  upstream  volume  stagnation  pressure 
and  Junction  stagnation  enthalpy  and  the  critical  mass  flux 
Is  returned. 


1 


r 


1 
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Assumptions  implicit  in  the  use  of  this  model  are: 

(1)  System  flow  is  isentroplc  and  one ~d imens  lonal  , 

(2)  The  liquid  phase  is  Incompressible, 

(3)  The  transfer  of  heat  and  mass  between  phases  is 
negligible  , 

(4)  Slip  between  phases  is  negligible  and, 

(5)  At  a Junction  the  vapor  phase  expansion  can  be 
described  by  a polytropic  process. 

The  Henry-Faus  k.e  model  is  applicable  in  the  subcooled  , 
saturated,  and  superheated  regions  with  stagnation  pressure 
limits  of  1 PSIA  to  2400  PSIA. 

2.5  THE  HOMOGENEOUS  EQUILIBRIUM  MODEL  (HEM)  (13) 

This  is  essentially  the  same  as  the  Sonic  model 
described  in  section  2.2,  page  15.  It  differs  only  in  how 
it  is  utilized  by  the  RELAP4  code.  The  HEM  critical 
velocity  and  the  Sonic  model  critical  velocity  are 
calculated  using  the  same  formula,  equation  (6).  The  HFM 
critical  mass  flux  is  then  calculated  from, 


G - { 0«>j 


(44) 


The  difference  between 
density  is  calculated.  In 


the  two  models  lies  in  how  the 
the  Sonic  model  this  density  is 


computed  from  the  upstream  volume  density  by  adjusting  for 
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frictional  losses  and  kinetic  energy  changes  to  the  Junction 
and  by  isentropic  expansion  fron  the  Mach  number  just 
upstream  of  the  Junction  to  Mach  one  in  the  Junction.  In 
HEM  the  density  is  calculated  by  assuming  that  the  upstream 
volume  density  and  the  pressure  and  enthalpy  returned  from 
the  steam  tables  are  stagnation  properties.  The  junction 
density  used  in  equation  (44)  is  obtained  by  an  isentropic 
expansion  from  the  upstream  volume  assumed  stagnation 
conditions  to  sonic  velocity  in  the  Junction. 

The  assumptions  Inherent  in  the  HEM  are  the  same  as  for 
the  Sonic  model  given  on  p a g,  e 16  with  the  exceptions 
mentioned  in  the  previous  paragraph. 

The  HEM  critical  mass  flux  is  also  tabulated  for  use  in 
RELAP4  for  input  values  of  upstream  volume  stagnation 
pressures  and  junction  stagnation  enthalpy.  The  Homogeneous 
Equilibrium  model  is  applicable  from  1 PSIA  to  3000  PSIA. 

2.6  THE  MODIFIED  MOMENTUM /HOMOGENEOUS  EQUILIBRIUM  MODEL  (14) 

The  Modified  Momen t um /Horn o g en eous  Equilibrium  model 
(MM/HEM)  is  the  same  as  the  HEM  model  except  when  the 
quality  is  less  than  a user  specified  transition  quality,  X 

I 

(the  default  value  for  this  quality  is  X T ■ 0.02). 

As  mentioned  in  Chapter  One,  the  simple  momentun  or 
inertial  flow  solution  is  quite  accurate  at  low  Junction 
pressure  ratios  but  tends  to  exaggerate  the  flow  rate  as 
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critical  conditions  are  approached.  in  order  to  compensate 
for  this  exaggeration  the  modified  momentum  model  utilizes  a 
correction  to  the  downstream  pressure  based  upon  the 
retardation  of  vapor  formation  by  the  surface  tension  of  the 
flashing  fluid. 

This  more  realistic  downstream,  or  back  pressure,  is 
calculated  from  the  empirical  relationship. 


If,  however,  the  modified  back  pressure  is  less  than  or 
equal  to  the  downstream  volume  static  pressure  then  the  flow 
estimate  is  based  on  inertial  considerations  only.  In  the 
transition  region  (0.0£  X < XT ) the  critical  flow  estimate  is 
taken  as  the  minimum  of  the  flow  rate  predicted  by  equation 
(47)  and  the  flow  rate  calculated  from. 


WMM/HEM  “ GHEMA(XT/XUP) 

The  assumptions  and  limitations  are  the  same  as  those 
for  the  HEM  with  the  additional  empirical  assumption  of 
equation  (45). 

2.7  JUNCTION  FLOW  RATE  SOLUTION  LOGIC 

If  one  of  the  critical  flow  models  Is  selected  for  a 
particular  Junction  the  flow  rate  solution  process  for  that 
Junction  Is  accomplished  as  follows. 

An  inertial  Junction  flow  rate  estimate  Is  obtained 
from  the  following  linearized  expression  derived  from  the 
Integrated  form  of  the  momentum  equation; 

Vf  - W + g (iP  - MU2)  it/I  (44) 
t + t c t 

where  all  of  the  parameters  are  evaluated  at  the  Junction. 

This  estimate  is  then  compared  with  that  for  the 
particular  critical  flow  model  selected  and  the  minimum 
identified.  The  fluid  acceleration  and  Junction  friction 
required  to  force  the  actual  calculated  flow  rate  to  the 
above  Identified  minimum  estimate  is  determined  and  the 
appropriate  terms  in  the  momentum  equation  are  nodlfled 
before  simultaneous  solution  with  the  continuity  and  energy 
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(48) 


equations. 


1 

3 I 

It  la  Important  to  note  that  regardless  of  which 
critical  flow  model  Is  selected,  a comparison  with  inertial 
flow  Is  always  effected,  regardless  of  the  conditions  In  the 
Junction  and  the  upstream  volume.  Thus,  If  It  is  known 
apriori  that  a particular  Junction  will  not  attain  critical 
conditions,  then  selecting  a critical  flow  model  for  that 
Junction  serves  only  to  increase  computer  processing  time. 
Additionally  if,  for  example,  the  Moody  critical  flow  model 
is  selected  for  one  of  the  Junctions  connecting  core 
volumes,  when  the  fluid  in  those  core  volumes  becomes 
superheated  near  the  end  of  the  blowdown  phase,  the  problem 
will  abort  because  the  Moody  model  is  not  applicable  in  the 
superheated  regime. 


r 
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3.  OPTIMUM  AREA  RATIO  DETERMINATION 

3.1  INTRODUCTION 

As  reported  in  Section  1.5,  the  RELAP4  code  in  its 
present  configuration,  incorrectly  computes  stagnation 
properties.  In  order  to  alleviate  this  problem,  the 
upstream  volume  flow  area  must  be  increased  in  order  to 
prevent  supersonic  velocities  from  occurring  in  the  volume 
Immediately  upstream  from  the  break  junction. 

This  chapter  presents  the  results  of  a study  done  to 
determine  the  optimum  area  ratio  (AR)  (ie.  break  junction 
flow  area/  upstream  volume  flow  area)  for  this  problem.  The 
optimum  area  ratio  was  found  to  be  D.7  (ie.  the  flow  area 
for  volumes  V5  and  V 2 1 must  be  increased  from  3 . A 7 to  5 . 2 & 
square  feet).  In  the  chapters  that  follow,  junctions, 
volumes,  and  heat  slabs  will  be  referred  to  hy  the  letters 
J,  V,  and  S,  respectively,  followed  by  the  number  of  that 
Junction,  volume,  or  heat  slab  as  defined  in  Figure  1 page 
7.  For  example,  J27  refers  to  the  break  junction  on  the 
vessel  side  of  the  break,  V5  refers  to  the  volume 
immediately  upstream  of  break  Junction  J27,  and  V21  refers 
to  the  volume  immediately  upstream  of  break  junction  J2R. 


3.2  GENERAL 


This  study  made  use  of  the  Hope  Greek  BWR  system  as 
modeled  for  RELAP4  by  EG&G  Idaho,  Inc.  with  Evaluation  Model 
controls  applied  (RELAP4-EM).  The  input  deck  for  this 
Evaluation  Model  is  presented  in  Appendix  A.  The  critical 
flow  model  required  for  the  Evaluation  Model  is  the  Henry- 
Fauske  model  in  the  subcooled  reg  line  and  the  Moody  model  in 
the  saturated  regime.  The  other  Evaluation  Model  controls 
are  irrelevant  to  this  area  ratio  study  and  discussion  of 
these  controls  will  be  delayed  until  Chapter  Four. 

In  order  to  understand  the  system  transient  response  to 
be  described  in  the  following  section,  the  control  action 
inputs  used  in  the  problem  must  be  explained. 

The  problem  starts  at  time  0.001  seconds  when  the  brent' 
occurs.  At  this  time  the  break  Junctions  J27  (vessel  side) 
and  J28  (pump  side)  are  opened.  At  time  0.002  seconds,  .12  5 
(the  pre-break  recirculation  flow  Junction)  closes  to 
complete  the  200S  break.  Also  at  this  time  the  feedwater 
valve,  J30,  closes,  the  two  recirculation  pumps,  modeled  by 
volumes  V6  and  V9,  are  shut  off  and  begin  to  coast  down,  and 
a reactor  SCRAM  is  initiated.  At  one  second  into  the 
transient  the  main  steam  stop  valve,  J31,  begins  to  close. 
This  is  a slow  closing  valve  and  is  not  fully  closed  until 
four  seconds  into  the  transient.  The  core  spray  system 
(CSS),  J32,  and  the  low  pressure  coolant  injection  system 


(LPCIS),  J33,  actuate  26.5  seconds  after  the 


mixture  level 


in  the  lower  downcomer,  V4,  reaches  21.5  feet.  As  a 
reference,  the  bottom  of  the  lower  downcomer  is  7.08  feet 
below  the  bottom  of  the  active  core.  Flow  through  these 
systems  does  not  begin,  however,  until  the  volumes  to  which 
they  are  discharging  (VI  for  the  CSS  and  V10  for  the  LPCIS) 
depressurize  to  304  PSIA  and  310  PSIA,  the  rated  discharge 
pressures  of  the  CSS  and  LPCIS  pumps,  respectively.  The 
automatic  depressurization  system  (ADS),  J26,  is  not 
activated  until  120  seconds  after  the  mixture  level  in  the 
lower  downcomer  reaches  21.38  feet  and  thus  is  not  a fc.ctor 
in  this  study. 

In  order  to  determine  the  optimum  area  ratio,  six 
RELAP4  runs  were  completed  with  area  ratios  of  1.0,  0.9, 
0.8,  0.7,  0.6,  and  0.5.  The  temporal  response  of  seven  key 

parameters,  listed  below,  to  the  Imposed  LOCA  transient  was 
studied; 

1)  Break  flow  rate,  vessel  side  of  break-J27, 

2)  Break  flow  rate,  pump  side  of  break-J28, 

3)  Core  inlet  flow-J29, 

4)  Lower  d own c ome r- V4 , mixture  level, 

5)  Lower  plenun-Vll,  pressure, 

6)  Clad  surface  temperature  in  the  hottest  (center) 
heat  slab  of  the  core-S23,  and 

7)  Heat  transfer  coefficient  of  the  fluid  adjacent  to 
the  hot  slab  of  the  core-S23. 
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In  general,  the  temporal  parametric  responses  were 
bounded  by  the  1.0  and  the  0.7  area  ratio  runs.  For 
purposes  of  clarity  only  these  two  cases  will  he  compared 
for  all  of  the  seven  parameters  considered.  For  break  flow, 
J27,  and  core  inlet  flow,  J29,  all  six  cases  will  be 
compared  for  the  time  Interval  where  the  greatest  (and  most 
Important)  differences  occurred. 

In  all  of  the  figures  presented  In  this  chapter,  the 
curves  are  Identified  by  numbers  (1.0,  0.9,  0.8,  etc.)  which 
correspond  to  the  area  ratio  used  to  produce  the  data  for 
that  curve.  The  data  used  to  produce  the  figures  In  thi3, 
and  subsequent  chapters,  was  extracted  from  the  PELAP4 
Plot/Restart  tape  rather  than  from  the  printed  output.  A 
Fortran  IV  program  written  to  extract  this  data  from  the 
tape  Is  presented  In  Appendix  B. 

3.3  ANALYSIS  OF  RESULTS 

The  most  Important  parameter  with  respect  to  blowing 
down  of  the  reactor  vessel  Is  the  rate  at  which  fluid  exits 
the  break.  For  a 200T  break  there  are  two  break  flows  of 
significance;  the  flow  from  the  vessel  side  of  the  break 
(J27)  and  the  flow  from  the  pump  side  of  the  break  (J28). 

Figure  3 shows  the  flow  rate  out  of  the  vessel  side  of 
the  break.  When  the  break  opens,  the  flow  rate  Increases 
sharply,  then  decreases  when  the  flow  conditions  become 


critical,  or  choked 
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As  the  escaping  fluid  Increases  In 
steam  quality,  the  flow  rate  decreases  slightly  until  the 
main  steam  line  closes  at  four  seconds.  The  flow  rate  Chen 
increases  momentarily.  The  flow  rate  remains  relatively 
constant  until  the  lower  downcomer  empties,  allowing  steam 
to  escape  from  the  break.  As  can  be  seen  from  Figure  4,  the 
lower  downcomer  empties  at  about  eleven  seconds  for  the  AR  - 
0.7  case  and  at  about  14.5  seconds  for  the  AR  - 1.0  case. 
These  times  correspond  to  the  sharp  decreases  in  break  flow 
seen  in  Figure  3.  Following  the  emptying  of  the  lower 
downcomer,  the  flow  rate  decreases  steadily  due  to  the 
gradual  equlllzatlon  of  pressure  between  the  containment  and 
the  reactor  vessel. 

The  inset  to  Figure  3 reveals  that  the  AR  « 1.0  case 
initially  peaks  at  a higher  flow  rate  than  does  the  AR  » 0.7 
case.  As  will  he  seen,  similar  results  were  obtained  for 
initial  break  flow  on  the  pump  side  of  the  break  (Figure  6, 
page  42).  A satisfactory  explanation  for  this  peculiar 
result  could  not  be  found.  While  the  difference  in  this 
initial  break  flow  is  substantial,  it  occurs  for  such  a 
short  time  that  it  is  insignificant. 

As  can  be  seen,  the  major  differences  between  the  two 
runs  lie  in  the  value  of  the  choked  flow  rate  attained 
immediately  following  break  initiation.  ror  essentially  the 
same  driving  pressure  differential  the  AR  - 0.7  case  yields 


a break  flow  rate  of  about  23,000  LBM/SEO  compared  to  only 
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16,000  LBM/SEC  for  the  AR  ■ 1.0  case.  This  results  in  the 
lower  downcomer  emptying  faster  for  the  AR  - 0.7  case. 
After  the  lower  downcomer  empties,  the  two  flow  rates  remain 
essentially  the  same,  with  the  AR  ■ 1.0  case  being  slightly 
greater  due  to  the  higher  Junction  pressure  differential 
which  in  turn  is  due  to  a lower  break  flow  rate  in  the 
initial  fourteen  seconds  of  the  blowdown. 

Figure  5 shows  the  lower  plenum  pressure  which  is 
essentially  the  driving  pressure  for  the  break  flow.  This 
pressure  decreases  slowly  during  the  first  four  seconds  of 
the  transient  because  the  loss  of  water  from  the  break  and 
the  loss  of  steam  out  of  the  main  steam  line  are  Just 
slightly  greater  than  can  be  compensated  for  by  the 


production  of  steam  in  the  core.  After  the  main  steam  line 
closes  at  four  seconds  the  pressure  increases  slowly  until 
emptying  of  the  lower  downcomer  allows  steam  to  escape  from 
the  break.  From  here  on  the  pressure  decreases  slowly  as 
the  reactor  vessel  and  containment  pressures  equalize. 

The  differences  between  the  two  cases  are  due  mainly  to 
the  time  it  takes  for  the  lower  downcomer  to  empty.  Thus 
the  AR  * 1.0  case  yields  a higher  maximum  pressure  due  to 
the  Increased  time  between  main  steam  line  shutoff  and 
emptying  of  the  lower  downcomer.  Following  the  initiation 
of  steam  flow  out  of  the  break  the  two  cases  show 


depressurization  rates  which  are  essentially  the  same. 


TIMF. 
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The  flow  out  of  the  pump  side  of  the  break  (J27)  is 
shown  In  Figure  6.  It  should  be  noted  that  J28  is  modeled 
for  a positive  flow  Into  the  system.  For  convenience. 
Figure  6 is  plotted  showing  positive  flow  out  of  the  break. 
The  flow  characteristics  on  the  pump  side  of  the  break  are 
significantly  different  from  those  on  the  vessel  side  of  the 
break  due  to  the  presence  of  the  recirculation  pump  which 
begins  to  shut  down  at  0.002  seconds.  The  break  flow 
initially  increases  at  the  same  rapid  rate  as  did  the  vessel 
side  of  the  break  and  then  decreases  rapidly.  This  is 
because  the  flow  is  being  restricted  by  the  recirculation 
pump  and  the  Jet  pump  drive  nozzles.  After  about  1.5 
seconds  the  break  flow  (J28)  decreases  slowly  and  is 
oblivious  to  other  changes  in  the  reactor  system. 

As  can  be  seen,  the  differences  between  the  two  cases 
shown  in  Figure  6 are  slight.  The  inset  in  Figure  6 shows 
that  the  peak  flow  rate  is  about  20Z  higher  for  the  AR  - 1.0 
case.  As  pointed  out  earlier,  this  difference  is 
insignificant  because  it  occurs  for  such  a short  period  of 
time. 

In  all,  it  is  apparent  that  the  area  ratio  change  had 
very  little  effect  on  the  flow  rate  from  the  pump  side  of 
the  break.  This  is  due  to  the  lower  break  junction  velocity 
caused  by  the  flow  restrictions  of  the  r e c i r c u 1 a t i on  pump 
and  jet  pump  drive  nozzles. 


FIGURE  6 

BREAK  FLOW  PUMP  SIDE  OF  BREAK  (J28) 
EVALUATION  MODEL 


Figure  7 shows  the  core  inlet  flow  (J29).  Immediately 
after  the  break,  the  core  flow  begins  to  decrease  as  the 
recirculation  pump  on  the  side  opposite  the  break  coasts 
down  and,  as  fluid  normally  flowing  vertically  upward 
through  the  core,  is  diverted  to  the  break  through  the  jet 
pump  on  the  break  side*  This  steady  decrease  continues 
until  the  Jet  pump  intakes  are  uncovered.  The  Jet  pump 
intakes  are  located  at  the  tops  of  volumes  V19  and  V20,  and 
uncover  when  the  mixture  level  in  the  lower  downcomer 
reaches  15.7  feet.  Reference  to  Figure  4 , page  3 R , reveals 
that  the  Jet  pump  Intakes  uncover  at  about  eight  seconds  for 
the  AR  » 0.7  case  and  at  about  ten  seconds  tor  the  AR  « 1.0 
case.  These  times  correspond  to  the  sharp  decreases  in  core 
flow  seen  in  Figure  7.  Uncovering  the  Jet  pump  intakes 
causes  the  remaining  pressure  differential  across  the  core 
(due  to  the  recirculation  pump  in  the  intact  loop  which  is 
still  coasting  down)  to  diminish,  and  the  core  flow  then 
reverses.  Ac  the  time  of  core  flow  reversal,  any  liquid 
remaining  in  the  core  is  replaced  by  steam  from  the  upper 
plenum  and  the  core  dries  out.  The  core  flow  remains 
negative  until  the  water  in  the  lower  plenum  begins  Co  flash 
to  steam  as  the  system  pressure  drops  to  below  saturation. 

Flashing  of  the  water  in  the  lower  plenum  causes  the 
lower  plenum  pressure  to  decrease  at  a slower  rate  than  does 
the  upper  plenum  pressure,  thus  the  core  inlet  flow  reverses 
again.  The  core  flow  then  decreases  as  the  liquid  Inventory 
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in  the  reactor  system  decreases.  The  flow  again  turns 
negative  and  remains  near  zero  for  the  remainder  of  the 
transient • 

The  differences  between  the  two  cases  can  again  he 
attributed  to  the  differences  in  the  initial  break  flow 
rates  discussed  earlier.  The  lower  break  flow  rate  for  the 
AR  » 1.0  case  causes  the  Jet  pump  drive  nozzles  to  uncover 
later  (ten  seconds  versus  eight  seconds),  and  thus  the  first 
core  flow  reversal  occurs  later.  Lower  plenum  flashing 
starts  when  the  system  pressure  drops  below  about  1000  PSIA, 
the  saturation  pressure  for  about  545  degrees  F (the  average 
water  temperature  in  the  lower  plenum)  . Referring  to  Figure 
5,  page  40,  the  system  pressure  drops  helow  1000  PSIA  at 
about  11.5  seconds  for  the  AR  ■»  0.7  case  and  at  about 
sixteen  seconds  for  the  AR  ■ 1.0  case.  These  times 
correspond  to  the  second  core  reversals  seen  in  Figure  7. 

The  temporal  response  of  clad  surface  temperature  is 
shown  in  Figure  8.  The  data  for  this  graph  is  obtained  by 
taking  the  surface  temperature  of  the  hottest  heat  slab 
(S23),  which  is  in  the  center  of  the  core.  Tt  should  he 
noted  that  this  system  is  modeled  with  the  core  heat  slabs 
representing  an  average  fuel  assembly.  Thus  the  surface 
temperature  of  heat  slab  S23  represents  the  hottest  fuel 
clad  temperature  of  an  average  fuel  assembly. 


A 

no  r e 
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core 

model  , 

such 

as  the  model  hv 

Hendrix 

(9)  , 

divides 
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core 
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two  regions,  one 
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representing 

the  average 

fuel  assembly,  and 

the 

other 

representing 

the  hot 

fuel  assembly.  While 

this 

more 

sophisticated 

model 

will 

more  accurately  predict 

core 

tempe  natures  , 

It  will 

also 

Increase  computer  processing 

time 

and  would  have  added  little  to  this  analysis,  which  Is 
concerned  with  relative,  rather  than  absolute,  parametric 
responses • 

Figure  8 Indicates  that  Initially  the  clad  surface 
temperature  follows  the  fluid  saturation  temperature 
(approximately  565  degrees  F ) closely  until  the  steam 
quality  In  the  core  becomes  appreciable.  At  the  time  of 
first  core  flow  reversal  the  temperature  of  the  cladding 
surface  Increases  sharply  as  the  core  dries  out  and  the  heat 
transfer  coefficient  of  the  fluid  adjacent  to  the  heat  slah 
decreases.  Figure  9 shows  the  heat  transfer  coefficient 
(HTC)  of  the  fluid  adjacent  to  S23.  The  HTC  decreases 
slowly  until  the  first  core  flow  reversal  where  It  drops  to 
nearly  zero  due  to  dryout  of  the  core.  The  time  for  this 
sharp  decrease  In  the  HTC  corresponds  to  the  time  of  first 
core  flow  reversal  (Figure  7)  and  to  the  sharp  Increase  In 
clad  surface  temperature  (Figure  8).  Fefering  hack  to 
Figure  8,  the  clad  surface  temperature  continues  Its  rapid 
increase  until  the  time  of  lower  plenum  flashing  and  second 
core  flow  reversal.  There  It  levels  out  somewhat  due  to  the 
liquid  flow  through  the  core  caused  by  lower  plenum 
flashing.  Following  the  flow  spike  caused  by  lower  plenum 
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flashing  the  surface  temperature  of  the  cladding  again 
increases,  levels  out,  and  remains  relatively  constant  until 
the  core  spray  system  (CSS)  becomes  effective.  (The  CSS 
does  not  become  effective  in  the  forty  seconds  of  transient 
time  shown ) . 

As  shown  by  Lalra  (19)  and  Bruch  (20)  if  the  transient 
is  allowed  to  progress  beyond  forty  seconds,  the  clad 
surface  temperature  will  gradually  decrease  until  ahout  l 2 S 
seconds,  when  it  will  begin  to  increase.  The  peak  cladding 
temperatures  will  not  be  reached  until  well  into  the  reflood 
phase  of  the  accident.  RELAP4-FL00D  is  not  available  for 
BWR  s and  therefore  the  peak  cladding  temperature  cannot  he 
predicted  using  RELAP4  in  its  present  configuration. 

Again,  the  differences  between  the  A R ■ 0.7  clad 
surface  temperature  case  and  the  AR  ■ 1.0  case  lie  in  the 
initial  break  flow  differences  which  directly  influence  the 
time  of  first  core  flow  reversal  (Figure  7).  As  pointed 
out,  the  times  of  first  core  flow  reversal  correspond  to 
drying  out  of  the  core  and  to  the  resultant  sharp  Increases 
In  cladding  surface  temperature  seen  in  Figure  R.  Following 
first  core  flow  reversal  the  clad  surface  temperature  for 
both  cases  follow  nearly  the  same  trends  with  the  AR  « 1.0 


case  remaining  below  the  AR 


0.7  case. 
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3.4  OPTIMUM  AREA  RATIO  SELECTION 

Figure  10  shows  the  break  flow  f ron  the  vessel  side  nf 
the  break  for  each  of  the  area  ratio  cases  studied.  In 
order  to  accentuate  the  differences  between  the  cases,  the 
transient  response  is  only  shown  for  the  transient  time  f ron 
one  to  ten  seconds. 

As  can  be  seen  from  Figure  10,  increasing  the  upstream 
volume  flow  area  (ie.  decreasing  the  area  ratio)  results  in 
increasing  the  break  flow  rate,  with  the  slight  exception  of 
the  time  period  from  five  to  ten  seconds  for  the  AR  « 0.5, 
0.6,  and  0.7  cases.  While  these  three  cases  are  different 
for  most  of  the  transient  shown,  they  never  differ  by  more 
than  four  percent,  and  for  all  practical  purposes  can  he 
considered  to  be  the  same. 

Similar  results  are  obtained  for  core  inlet  flow  shown 
in  Figure  11  for  transient  times  from  two  to  eleven  seconds. 
The  rain  point  of  interest  is  the  time  of  first  core  flow 
reversal  which  occurs  earlier  as  the  area  ratio  is  decreased 
until  the  area  ratio  reaches  0.7.  Further  reductions  in 
area  ratio  produce  no  significant  changes. 

As  pointed  out  in  the  introduction  to  this  chapter,  the 
parametric  responses  studied  generally  fell  between  the  AR  » 
1.0  case  and  AR  • 0.7  case.  It  was  also  noted  that  results 
similar  to  those  demonstrated  by  Figures  10  and  11  were 
observed  for  all  of  the  parameters  studied  (le.  decreasing 
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the  area  ratio  below  0.7  resulted  in  no  significant  change 
In  the  response  of  the  parameter).  For  these  reasons,  an 
area  ratio  of  0.7  was  chosen  as  the  optimum  value  for  this 
problem. 

It  is  important  to  note  that  using  the  RF.LAP4  code 
without  correcting  for  the  stagnation  properties  calculation 
error  (ie.  using  AR  - 1.0)  yields  results  which  are  not 
conservative.  That  is,  low  break  flow  rates  and  low  clad 
surface  temperatures  are  predicted. 
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4.  CRITICAL  FLOW  MODEL  SENSITIVITY  STUDY 

4.1  I NTRODUCTION 

A sensitivity  study  of  the  critical  flow  models 
discussed  in  Chapter  Two  using  the  Standard  Model  (best 
estimate)  configuration  of  the  RELAP4  code  (RELAP4-SM)  was 
made.  A comparison  was  also  made  with  the  Evaluation  Model 
(RELAP4-EM).  It  should  be  recalled  that  the  critical  flow 
model  for  the  Evaluation  Model  is  the  Henry-Fauske  model  in 
the  subcooled  regime  and  the  Moody  model  in  the  saturated 
regime.  The  input  data  decks  used  in  this  study  are 
discussed  in  Appendix  A. 

The  results  of  the  area  ratio  study  presented  in  the 
previous  chapter  were  applied  to  the  sensitivity  study  runs. 
That  is,  the  flow  areas  for  volumes  V5  and  V 2 1 (the  volumes 
immediately  upstream  from  the  break  Junctions,  J 2 7 and  J 2 ■*  , 
respectively)  were  increased  from  3.57  to  5.24  square  feet. 
The  Evaluation  Model  run  presented  in  this  chapter  is  the 
-0.7  case  from  Chapter  Three. 

4.2  ANALYSIS  OF  RESULTS 

Figure  12  shows  the  break  flow  rate,  for  the  first 
twenty  seconds  of  the  transient,  from  the  vessel  side  of  the 
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break  (J27)  for  the  six  cases  studied-  The  transient 
response  for  this  parameter  is  uninteresting  from  twenty  to 
forty  seconds  and  is  similar  to  the  response  shown  in  Figure 
3,  page  36.  In  Figure  12,  and  in  the  figures  and  discussion 
that  follows,  the  different  cases  are  Identified  by  the 
symbology  presented  in  Table  2. 

TABLE  2.  CRITICAL  FLOW  MODEL  SYMBOLOGY 


CRITICAL  FLOW  MODEL 

SYMBOL 

SONIC 

S 

HOMOCENEOUS  EQUILIBRIUM 

HFM 

MODIFIED  MOMENTUM/HOMOGENEOUS  EQUILIBRIUM 

MMHF.V 

HENRY-FA USKE 

HF 

MOODY 

M 

EVALUATION  MODEL 

FM 

Figure  12  reveals  that  the  EM  yields  the  highest 
average  break  flow  rate  for  the  initial  period  of  the 
blowdown,  and  is  followed  by  the  Moody,  the  He n r y -Fa  us  V e , 
the  MMHEM,  the  HEM,  and  the  Sonic  models. 

With  the  execption  of  the  Moody  model,  the  results 
depicted  in  Figure  12  are  as  would  have  been  expected  based 
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on  the  discussion  on  break  flow  given  in  Chapter  Three. 
That  is,  following  the  rapid  decrease  in  break  flow 
characteristic  of  all  of  the  models  at  about  eleven  seconds, 
the  Moody  model  exhibits  a peculiar  Jump  in  hreak  flow  at 
about  11.5  seconds. 

Analysis  of  the  Moody  data  indicates  that  the  lower 
downcomer  empties  at  11.2  seconds,  allowing  steam  to  escape 
from  the  break.  At  11.6  seconds,  the  lower  downcomer 
refills  to  a level  of  about  1.2  feet,  and  the  flow  out  of 
the  break  again  becomes  mostly  liquid.  The  lower  downcomer 
empties  again  at  12.6  seconds,  and  remains  empty  for  the 
remainder  of  the  transient.  This  oscillation  of  lower 
downcomer  water  level  is  caused  by  increased  steam 
generation  in  the  core  starting  at  about  11.4  seconds.  The 
increased  steam  generation  forces  water  from  the  lower 
plenum,  through  the  Jet  pump  intakes,  and  Into  the  lower 
downcomer.  The  increase  in  core  steam  generation  rate  is 
caused  by  a return  to  nucleate  boiling  following  first  core 
flow  reversal  and  core  dryout.  The  return  to  nucleate 
boiling  will  he  explained  in  the  discussion  of  the  heat 
transfer  coefficient  which  follows. 

The  oscillation  of  lower  downcomer  mixture  level, 
described  above  for  the  Moody  model,  does  not  occur  for  the 
other  five  cases. 

The  hreak  flow  from  the  pump  side  of  the  break  ( J 2 0 ) 
shoved  only  slight  differences  (less  than  two  percent) 


between  the  six  cases  studied 
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All  six  cases  were  very 
similar  to  the  response  shown  In  Figure  6 , pap,  e 42,  and  thus 
are  not  presented. 

Table  3 shows  the  times  for  uncovering  the  jet  punp 
drive  nozzles  and  the  times  for  emptying  of  the  lower 
downcomer  for  the  six  cases  studied.  (Recall  from  Chapter 
Three  that  the  jet  pump  drive  nozzles  uncover  when  the 
mixture  level  In  the  lower  downcomer  reaches  1 5 . 7 feet, 
which  corresponds  to  a mixture  level  of  about  8.8  feet  In 
the  active  core).  The  times  presented  In  Table  3 are  a 
direct  result  of  the  initial  break  flow  rates  depicted  In 
Figure  12. 


TARLF. 

3.  TIMF.  TO 

UNCOVER  JET  PUMP  DRIVE 

NOZZLFS  AND 

TO 

EMPTY  LOWER  DOWNCOMER 

CRITICAL 

MODEL 

FLOW  TIMF(SEC)  TO  UNCOVER  TIME(SEC)  TO  FM  PT  V 

JET  PUMP  DRIVE  NOZZLFS  LOWFR  DOWNCOMFR 

EM 

7.  5 

1 0.  8 

M 

7.  6 

11.2 

HF 

7.  7 

11.4 

MMHEM 

8.  1 

1 2 . a 

HEM 

8 . 8 

14.0 

S 


9.  S 


1 4.  8 


J 


r 
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As  discussed  In  Chapter  Three,  the  time  that  the  Jet 

i 

punp  drive  nozzles  uncover  should  correspond  to  the  time  of 
first  core  flow  reversal,  and  the  time  that  the  lower 
downcomer  empties  should  correspond  to  Che  time  of  second 
core  flow  reversal.  This  is  indeed  Che  case,  as  can  he  seen 
in  Figure  13,  which  shows  the  core  Inlet  flow  (J20)  for  the 
six  cases  studied  from  zero  to  twenty  seconds  into  the 
transient.  The  remainder  of  the  transient  was  similar  to 
the  EM  response  shown  in  Figure  7,  page  44.  That  is,  the 
core  flow  remained  near  zero  with  small  oscillations  due  to 
varying  core  steam  generation  rates. 

The  transient  responses  seen  in  Figure  13  are  as  would 
have  been  expected  from  the  discussion  on  core  inlet  flow 
given  in  Chapter  Three,  again  with  the  exception  of  the 
Moody  model.  The  relatively  small  flow  spike  at  the  time  of 
second  core  flow  reversal  for  the  Moody  case  is  due  mainly 
to  the  return  to  nucleate  boiling  discussed  earlier.  This 
return  to  nucleate  boiling  results  in  increased  amounts  of 
steam  being  generated  in  the  core  which  counteracts  the 
effects  of  lower  plenum  flashing  on  core  differential 
pressure.  The  end  result  is  the  smaller  flow  spike  seen  at 
the  time  of  second  core  flow  reversal  in  Figure  13. 

Figure  14  shows  the  heat  transfer  coefficient  (UTC)  at 
the  hottest  heat  slab  in  the  core  for  the  six  cases  studied. 
It  should  he  noted  that,  while  the  HTCs  depicted  in  Figure 
14  are  for  the  hottest  heat  slab  in  the  core,  for  Che  HFM, 


Sonic,  and  Moody  models  the  hottest  heat  slab  turned  out  to 
be  heat  slab  S24  (the  slab  immediately  above  the  center  heat 
slab,  S23)  and  for  the  other  three  models  S23  was  the 
hottest  heat  slab.  This  result  is  explained  as  follows. 

One  of  the  main  differences  between  the  Fvaluatlon 
Model  (RELAP4-EM)  and  the  Standard  Model  (RE1AP4-SM)  is  In 
the  heat  transfer  logic  used  in  the  two  cases.  For  the  case 
of  RELAP4-FM.  it  is  stated  in  Appendix  K to  10  CFR  50  that 
"after  critical  heat  flux  is  first  predicted  at  an  axial 
fuel  rod  location  during  blowdown,  the  calculation  shall  not 
use  nucleate  boiling  heat  transfer  correlations  at  that 
location  subsequently  during  the  blowdown  even  if  the 
calculated  local  fluid  and  surface  conditions  would 
apparently  Justify  the  re-establishnent  of  nucleate 
boiling".  RELAP4-SM,  on  the  other  hand,  allows  re- 
establishment of  nucleate  boiling  if  the  conditions  warrant. 
For  the  HEM,  Sonic,  and  Moody  model  cases,  at  the  time  of 
first  core  flow  reversal,  the  conditions  were  such  that 
nucleate  boiling  was  re-established  at  heat  slab  S23,  while 
film  boiling  continued  at  beat  slab  S24.  Vucleate  boiling 
conditions  persisted  for  about  ten  seconds  at  heat  slab  S 2 3 
for  both  the  HEM  and  Sonic  models  and  for  about  three 


seconds  for 

the 

Moody  model.  This 

resul ted 

1 n 

the 

temperature  of 

heat 

slab  S23  being  lower 

than  that 

for 

heat 

slab  S24  in 

a 1 1 

three  cases  for  the 

remainder 

o f 

the 

transient. 


Figure  15  shows  the  cladding  surface  temperature  at  the 
hottest  heat  slab  In  the  core  for  the  six  cases  studied. 
Again,  for  the  HEM,  Sonic,  and  Moody  model  cases,  the 
hottest  heat  slab  was  S24,  while  for  the  other  three  cases 
S23  was  the  hottest. 

Table  4 shows  the  peak  cladding  surface  temperature 
attained  during  the  blowdown  for  each  of  the  six  cases 
studied . 


TABLE  4.  PEAK  CLADDING  SURFACE  TEMPERATURF 


CRITICAL  FLOW 
MODEL 


PEAK  CLADDING  SURFACE 
TEMPERATURE  (DEC  F) 


EM 

1340 

HF 

o 1 0 

009 

MMHEM 

RQ8 

HEM 

8 5 0 

S 

8 4 6 

As  can  be  seen  from  Figure  15  the  surface  temperature 
decreases  rapidly  for  the  Moody  model  case  at  about  eleven 
seconds.  This  decrease  corresponds  to  the  r e -e s t ah  1 1 s h men t 
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of  nucleate  boiling  which,  for  this  nodel  only,  occurs  In 
heat  slab  S24,  as  well  as  In  heat  slab  S23.  This  return  to 
nucleate  boiling  at  heat  slab  S24  is  characterized  by  the 
slight  Increase  in  the  HTC  seen  at  about  eleven  seconds  In 
Figure  14.  With  the  exception  of  this  one  peculiarity  for 
the  Moody  model,  the  surface  temperature  transient  responses 
for  all  six  cases  follow  trends  that  are  consistent  with  the 
discussion  of  clad  surface  temperature  In  Chapter  Three. 


4.3  CONCLUSIONS 

From  the  analysis  of  the  previous  section  It  Is 
apparent  that  the  Evaluation  Model,  as  expected,  produced 
the  most  conservative  results.  As  pointed  out,  this  Is 
primarily  due  to  the  different  heat  transfer  logic  In  the 
PELAP4-EM  and  RELAP4-SM  configurations. 

Rased  on  initial  break  flow  rate  (Figure  12,  p a p e 5R) 
■ it  is  clear  that  the  Moody  model  is  most  conservative, 

followed  by  He n r v -F a u s k e , “MHEM,  HEM,  and  Sonic.  Note, 
however,  that  the  Moody  and  Henry-Fa  uske  models  differ  by 
only  about  five  percent. 

Rased  on  peak  cladding  surface  temperature  attained 
(Table  4,  page  63)  one  would  conclude  that  the  Henry-Fauske 
model  is  more  conservative  than  the  Moody  model,  followed  by 
the  other  models  in  the  same  relative  order  as  above.  Note 
again,  that  the  difference  between  the  *'oody  and  Henry- 
Fauske  models  is  slight  (only  about  one  percent). 
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For  these  reasons,  the  Moody  and  Henry-Fauske  models 
are  rated  as  equivalent  based  on  conservative  results, 
followed,  In  order  of  their  relative  degrees  of 
conservatism,  by  the  MMHEM,  the  HEM,  and  the  Sonic  critical 
flow  mode  1 . 


SUMMARY 


5 . 


5.  1 RESULTS 


This  study  utilizes  the  R FLA P4  /MOP  5 Computer  Code  to 
analyze  the  effects  of  critical  flow  modeling  on  the 
t he r raa 1 -h y d r a u 1 i c transient  response  of  a Ceneral  Fleetrlc 
boiling  water  reactor  to  a major  primary  coolant  line 
rupture.  Included  in  the  study  is  a presentation  of  the 
equations,  assumptions,  and  limitations  of  the  critical  flow 
models  available  for  use  in  R E L A P 4 . Additionallv,  an 
evaluation  of  a temporary  solution  to  a RELAP 4 coding  error 
associated  with  stagnation  properties  calculation  is 
presented  . 

The  results  of  this  investigation  indicate  that: 

(1)  A solution  to  the  stagnation  properties 
c a 1 c u 1 a 1 1 ona  i error  that  exists  in  the  RELAP4  code,  when 
applied  to  the  Evaluation  Model,  (RELAP4-EM)  provides  a 
conservative  evaluation,  relative  to  the  Standard  Model 
(RELAP4-SM),  regardless  of  which  of  the  five  critical  flow 
models  available  in  the  code  is  selected; 

(2)  Of  the  five  critical  flow  models  available  for  use 
in  RELAP4,  the  Moody  model  and  the  Henry-Fauske  model  are 
nearly  equivalent  in  their  relative  degrees  of  conservatism. 
Following  these  models,  in  order  of  relative  degree  of 
conservatism,  are  the  'odified  ’'omen  t un  / Romo  g e n eo  us 


L 


Equilibrium  model,  the  Homogeneous  Equilibrium  model,  and 
the  Sonic  model  ; 

(3)  In  order  to  best  alleviate  the  stagnation 
properties  coding  error,  the  flow  areas  for  the  volumes 
Immediately  upstream  from  the  recirculation  line  rupture 
should  be  Increased  such  that  the  area  ratio  ^break  flow 
area  / upstream  volume  flow  area)  is  0.7. 

All  of  the  critical  flow  models  in  RELAP4,  with  the 
exception  of  the  Sonic  model,  require  the  use  of  stagnation 
properties  for  estimation  of  flow  rate  under  critical  or 
choked  conditions.  The  RELAP4  code,  in  its  present  form, 
incorrectly  computes  stagnation  properties  when  the  flow 
rate  is  such  that  the  code  predicts  supersonic  velocities  in 
the  volume  immediately  upstream  of  the  junction  where 
critical  flow  conditions  exist.  The  results  of  this 
investigation  indicate  that  the  optimum  solution  to  this 
problem  is  to  arbitrarily  increase  the  upstream  volume  flow 
area  such  that  the  area  ratio  is  0.7. 

The  Hope  Creek  boiling  water  reactor,  included  as  a 
sample  problem  in  the  RELAP4  Users  Manual,  was  used  to 
analyze  the  effects  of  critical  flow  modeling  on  the 
predicted  blowdown  response  of  this  BWR  to  a double-ended 
guillotine  rupture  of  the  primary  coolant  recirculation 
line.  Applying  the  above  mentioned  stagnation  properties 
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which  of  the  five  critical  flow  models  available  in  the  code 
is  selected.  Applying  the  "fix"  to  the  Standard  Model  it 
was  found  that  the  Moody  and  Henry-Fauske  critical  flow 
models  provide  nearly  equivalent  evaluations  in  terms  of 
peak  cladding  temperatures  attained  during  the  blowdown. 
The  other  three  critical  flow  models,  in  order  of  their 
relative  degrees  of  conservatism,  were  found  to  be  the 
Modified  Momen tum/Homogeneous  Equilibrium  model,  the 
Homogeneous  Equilibrium  model,  and  the  Sonic  model. 

5.2  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Based  on  this  investigation  the  following  recommenda- 
tions for  future  research  are  presented: 

(1)  While  this  Investigation  rated  the  critical  flow 
models  by  their  relative  degree  of  conservatism,  comparison 
of  the  results  of  RE LA  PA  analysis  with  LOFT,  SEMISCALE,  or 
other  large  scale  experimental  data  would  provide  valuahl; 
insight  as  to  which  of  the  models  most  accurately  predict 
flow  rate  under  critical  or  choked  conditions. 

(2)  The  heat  transfer  logic  employed  by  the  RF.LAPA  code 
appears  to  have  a very  significant  effect  on  the  overall 
performance  of  the  code.  Analysis  of  the  effects  of 
changing  heat  transfer  and  critical  heat  flux  correlations 
would  be  invaluable  to  future  users  of  the  code. 
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APPENDIX  A 


INPUT  DATA  DECKS 

This  appendix  provides  a listing  of  the  Input  data 
decks  for  the  Evaluation  Model  (RELAP4-EM)  and  the  Standard 
Model  (RELAP4-SM)  used  in  this  investigation.  The  data 
consists  of  title  cards,  data  cards,  and  comment  cards.  The 
comment  cards  are  prefixed  with  an  asterisk  (*). 

Table  5 is  the  Input  deck  for  RELAP4-EM  that  was  used 
in  the  area  ratio  study  presented  in  Chapter  Three.  The 
only  two  cards  that  were  changed  to  produce  the  different 
area  ratio  runs  were  the  volume  data  cards  050051  and 
050211.  Card  050051  is  is  the  data  card  for  volune  V 5 , and 


card  050 

211  is  the 

data 

card  for 

volume  V 2 1 . 

Th  e 

underl Ined 

number 

on  these 

cards 

is  the 

volume  flow 

area.  This 

quantity 

was  varied 

f r on 

3.57  square  feet  for 

the 

area  ratio 

CAR ) - 1 

.0  case  to 

7 . 34 

square 

feet  for  the 

AS 

* 0 . 5 case. 

The  example  shown  in  Table  5 is  the  AP  * 0.7  case. 

Table  6 is  the  input  deck  for  RE1.AP4-SM  used  in  the 
sensitivity  study  of  Chapter  Four-  Note  that  the  volune 
flow  areas  for  V 5 and  V21  (cards  O50051  and  O502111  are  5.24 
square  feet,  the  optimum  value  found  in  Chapter  Three.  The 
only  three  cards  changed  in  this  deck  to  produce  the 
different  critical  flow  model  runs  were  the  title  card  (the 
f Irsi  card  in  Table  6),  and  the  Junction  data  cards  0 0 0 2 7 l 
and  080281,  for  break  junctions  J 2 7 and  J2R,  respectlvelv. 
The  underlined  quantities  on  these  cards  are  the  critical 
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flow  model 

selection 

ind  ices 

ICHOKE  and 

JCHOKE 

. JCHOKE  is 

the  first 

underlined 

quantity 

, ICHOKE 

is  the 

second  . 

The 

values  of 

JCHOKE  and 

ICKOKE 

determine 

wh  i c h 

critical 

f 1 ow 

model  is  to 

be  used  for  that 

junction. 

Table 

V,  pane 

I-I  4 7 

of  the  R E l,  A P 4 Users  Manual  Is  used  to  determine  the  values 
of  JCHOKE  and  ICHOKE  for  the  desired  critical  flow  model. 
In  the  example  shown,  for  junctions  J27  and  J2P,  the  values 
are  JCHOK*  - 2,  ICHOKE  - 1 which  indicate  that  the  Sonic 


critical 

flow 

model  is  to 

b e 

used 

for 

J unc  t ions 

J 2 7 

and  J 2 P . 

Note 

that 

all 

other  junct 

Ions 

.ire 

set 

up  to  use 

the 

inertial 

flow 

mode 

1 ( 1 e 

. JCHOKE  - 

-i . 

ICHOKE  - 

1 ) . 

TABLE  5.  EVALUATION  MODEL  INPUT  DATA  DECK 
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• 3 UR  EVALUATION  MODEL  2 001  BREAK  USING  7X7  FUEL 
*2345678901234567890123456789012345678901234567890 

• 1 000  l L DM  EDI  NTS  TRP  VOL  SUB  TDV  JUN  PMP  CXV  N L 1C  FLL  SLAV  COM  MAT  COR 

010001  -2  9 5 10  23  4 0 3 3 2 4 2 4 2 5 1 0 5 5 0 1 

010002  3388.  1.0 

010003  0111 

020000  AP  11  JU  27  JW  23  JU  29  JO  IS  ML  4 SR  23  CR  23  TO  27 
030010  1 100  5 0 .001  .00001  .1  2000. 

030020  3 3 30  .01  .000001  1. 

030030  20  10  2 0 .01  .000001  60. 

030040  20  25  l 0 .01  .000001  130. 

030050  20  25  1 0 .01  .00001  2000. 

030003  50 

040010  l t 0 0 60.0  0. 

040020  1-4  3 0 40.0  0. 

040030  2100  0.0  0. 

040040  3100  0.001  0. 

040050  4100  .002  0. 

040060  5 1 0 0 1 . 0 0. 

040070  6100  1.E6  0. 

040080  7-5  4 0 21.5  26.5 

040090  8100  .001  1.E6 

040100  9-5  4 0 21.38  120. 

• 5 X X X X VOLUME  DATA  CARDS 


03001  1 

2 

0 

1 0 3 9 . 

-1.  .091632 

1 54  8 . 

1 2 . 

2 5 

22.25 

0 

2 2 4 . 

0. 

2 9.92 

0 

05002  l 

■» 

0 

1023.28 

- 1 . .002  44  8 

6 70  3 . 

12  . 

2 5 

5.27 

0 

4 9 7 . 

0. 

4 1.75 

0 

050031 

l 

0 

1024.98 

- 1 . . 969 

3 7 7 7 . 

2 1 . 

2 5 

2 1.25 

0 

-97. 

0. 

51.25 

0 

05004  l 

3 

0 

1031.93 

5 32  . 3 -1  . 

3177. 

32  . 

23 

32.2  3 

0 

165. 

0. 

10.02 

0 

05005 l 

0 

0 

102  4 . 

532.3  -1 . 

12  s . 2 5 

2 6 . 

2 6 . 

0 

5.24 

0. 

-14. 

0 

050061 

0 

0 

1126.17 

533.25  -1 . 

6 3 . 

3. 

3. 

0 

4 . 0 

0. 

-14. 

0 

0500  7 1 

0 

0 

1226.17 

533.25  -1 . 

3 66  . 

40  . 

53 

40.53 

0 

3 . 2 

0. 

-14. 

0 

0 5 00"  1 

0 

0 

1 02  4 . 

5 12  . 10  -1  . 

136. 

2 6 . 

2 6. 

0 

3.67 

<J. 

-14. 

0 

0 5 0 0 <5  1 

0 

0 

1126.17 

533.25  -1 . 

63  . 

3. 

3 . 

0 

4 . 0 

0 . 

-14. 

0 

050101 

0 

0 

1226.17 

5)3.25  -1 . 

366  . 

40. 

53 

40.53 

0 

3 . 2 

0. 

-14. 

0 

050111 

1 

0 

1362.16 

5 32  . 5 3 - 1 . 

2 13  1.5 

1 7. 

20 

17.20 

0 

12  0. 

0. 

0. 

: 3 

050121 

> 

0 

1054.67 

549.337  -l . 

22  3 . 

2.85 

2 . 

95  0 

81 

.091 

. 04  7 ) 

1 7.  fl 

2 1 3 

0 5 0 1 )1 

2 

0 

1052.97 

- 1 . . 03  7 1 5 7 

22). 

2 . 75 

•8 

L • 

75  0 

8 1 

.09  1 

. 04  n 

20.67  14 

0 50 1 4 t 

0 

1051.52 

- 1 . .05056  7 

8 1.09 

I . 

1 . 

0 

31 

.09  1 

. 04  7 3 

2 3.4 

2 1 5 

050151 

0 

104  6.  37 

- 1 . .09  1 994 

2 2 3 . 

2 . 75 

•j 

75  0 

3 l 

. C 9J 

. 04  7 3 

2 4 4 

2 1 6 

050161 

0 

1047.65 

-l.  . 1 3 98  7 5 

2 2 3 . 

2 . 75 

75  0 

3 1 

. 09  I 

. 04  73 

2 7.1 

7 1 

250171 

i 

0 

1054.28 

532.53  -l. 

942  . 4 7 

1 8. 

12 

19.12 

0 

52.30  1 3 . 

0. 

13 

0 5 0 l 3 i 

0 

1045. 36 

5 32  . 6 7 - l . 

534.53 

12.35  l 

2.35 

0 

4 3.63 

5 0 • 

17.92 

l 

050191 

0 

0 

1)59.2 

5 32  . 5 3 - 1 . 

115. 

15.7 

5 . 7 

0 

19.89 

0. 

10.02 

0 

050201 

0 

0 

1)59.2 

532.53  -1 . 

115. 

15.7 

1 5 

■j 

. / 

0 

19.69 

0. 

10.02 

0 

050211 

0 

0 

1330.3 

5 32  . 30  - 1 . 

10.75 

*»  _ 0 

2 

. 2 

0 

5.24 

0 . 

10  . 

0 

3 502  2 l 

«• 

0 

14.7 

100.  0.6 

3 . 4 3 E 6 

2 50  . 

0 

. 0 

0 

1 000  . 

100  . 

-10. 

0 

050231 

0 

1055.48 

532.53  -t . 

20.5 

1.25 

l 

. 2 5 

0 

81.09 

0. 

16.77 

1 2 

•6000  LIOL’IO  LEVEL  CALCULATION  CARD 
060000  II  23  12  13  14  15  16  17  IS  19  20  1 

360001  0.0  0.0  0.0  0.0  0.0 

•60002  UALLIS  SLI?  CORRELATION  CARD 
•60002  VALSCl  VALSC2 
* 6 X X X 1 RUBBLE  DATA  CARDS 
• 6 X X X 1 ALPH  VBUS 
0600  1 I 1.0  0.0 

060021  0.0  3.0 

060031  0.8  3.0 

06004  1 0 . 1 . E 6 

•60001  SLIP  VELOCITY  CARD 
•SXXXY  JUNCTION  DATA  CARD 

•2345473901234567890l23456789O1234567990l23456'f89Ol234567390l234567?9012 


0800  1 l 

1 

2 

0 

Cl 

2 3 5 8 3 

. 1 ’l . I 52.17  0 i 

. 0 

. 0. 

1 

0 

3 

0 

0.0 

0.  11 

0 

l . 

0 0 

03002  1 

2 

3 

0 

3 

-0  74  . 

3135  312.  52.25 

0 . 

0. 

3 • 

l 

0 

3 

0 3. 

.00.  11 

0 

0.0  0 

380031 

4 

0 

0 

2 5 4 4 4 

. 16  5 . 42  . 0 0 . 

A 0 

. 0. 

1 

0 

3 

0 

0.0 

0.  11 

0 

0 . 

0 0 

090041 

4 

20 

0 0 

10041. 

5 J. 112  25.  30. 

7 

084 

l . 

l 7 

2 

0 

0 2 

D.O  0. 

1 1 

0 0.0 

08005  1 

4 

2 l 

0 

0 

4 750. 

3.6  7 1 1.0  0. 

0 

. 0. 

1 

0 

3 

0 

0.0 

0.  1 1 

0 

0. 

0 0 

03006  l 

5 

6 

-l 

0 

4 750  . 

3.67-13.  0. 

0 

. 0. 

3 

0 

3 

0 

0.0 

0.  11 

o 

o . 

0 

0 9 00  ?t 

6 

* 

/ 

l 

0 

4 750  . 

3.2  -13.  0 . 

0 

. 0. 

3 

0 

3 

0 

0.0 

0.  11 

0. 

o 

. 0 

380081 

/ 

20 

0 

0 

4 750  . 

.538  25.  47.  . 

2 3 

73  6 

. 8 

0 

0 

0 

2 0. 

.0.11 

n 

0 

. 0 

09000  t 

4 

3 

0 

0 

4 750  . 

3.67  11.0  0. 

0 

. 0. 

1 

0 

3 

0 

0.0 

0.  11 

0 

0. 

0 

080101 

.9 

q 

0 

4 750  . 

3.67  -13.  0. 

0 

. 0. 

0 

0 

3 

0 

0. 

0.  1 1 

0 

0 . 

0 

0801  1 l 

9 

i 5 

2 

0 

4 750  . 

3.2  -13.  0. 

0 

. 0. 

0 

0 

3 

0 

0 . 

0.  1 1 

0 . 

0 

090121 

13 

l ) 

3 

0 

4 750  . 

.538  25.  47.  , 

.2  3?)  6 

. 9 

0 

0 

0 

■9  D 

. 0.  11 

0 

1 

r\ 

030131 

4 

1 9 

0 0 

1004  1 . 

5 3.112  25.  30. 

1 , 

084 

l 7 

: 

0 

0 1 

o . o . 

1 1 

0 

.0.0 

080141 

2 3 

» 9 

0 

0 

2 66  3** 

. 81.09  17.92 

0. 

0. 

0. 

0 

o 

3 

0 0 « 

.0.11 

o 

l 

.o  o 

330151 

12 

1 3 

0 

0 

2 46  39 

. 81.09  20.67 

0 . 

0. 

0. 

0 

0 

3 

0 3, 

. 0 . 11 

0 

l 

.0  1 

030161 

1 3 

l 4 

0 

0 

2 66  39 

. 81.0»  23.42 

0 . 

0. 

0. 

0 

0 

3 

0 0, 

, 0.  11 

0 

l 

.0  0 

TABLE  5-  (coat.) 


7 A 


080 l 7 t 

l 4 

l 5 

0 0 26639 

. 91 

. 09 

24.42  0 . 0 

. 0. 

0 

0 

3 

0 

0.  0.  11  0 1.0 

0 

080181 

l 5 

l ft 

0 0 26639 

. 91 

.09 

27.17  0 . 0 

. 0. 

0 

0 

3 

0 

0 . 0.  11  0 l .0 

0 

080191 

1 6 

l 

0 0 26639 

37 

. 13 

29.92  0.  0 

. 0. 

0 

0 

3 

0 

0 . 0 . 1 l 0 1.0 

0 

080201 

l 1 

l 7 

0 0 2 944  . 

2. 

246 

17.0  0.0 

. o. 

1 

0 

3 

0 

1.0  0. 

1 1 n o.  o o 

0802  1 1 

1 7 

18 

0 0 2944. 

26 

.20 

18.02  0.  0 

. 0. 

0 

0 

3 

0 

o.  o.  ii  o o.o 

0 

080221 

l 8 

1 

0 0 2944. 

43 

.635 

30.02  0.  0 

. 0 , 

0 

0 

3 

0 

0 . 0.  11  0 1.0 

0 

080231 

l 9 

1 l 

0 0 14791 

.5  19 

.69 

10.12  0.  0 

. 0. 

1 

0 

3 

0 

0.  0.  11  0 0.0 

0 

080241 

20 

l 1 

0 0 1 4 791 

.5  19 

. 89 

10.12  0.  0 

. 0. 

l 

0 

3 

0 

0.  0.  11  0 0.0 

0 

080251 

2 l 

5 

0 l 4750. 

3. 

6 7 

11.0  0.0 

. 0. 

0 

0 

3 

0 

0.  0.  It  0 0.0 

0 

0802  6 l 

3 2 

2 

0 3 0 . 

.484 

55.0  1 . 1 

. 1 . 

0 

0 

0 

0 

o.  o.  u n o.o 

0 

0802  71 

: i 

2 j 

0 4 0. 

3. 

67 

n.  o.i 

. . 5 

0 

0 

2 

0 

0.  .9  11  0 0.0 

0 

080291 

22 

5 

0 4 0 . 

3. 

87 

11.  o.  . 

5 l . 

0 

0 

•* 

0 

0.  .9  It  0 0.0 

0 

090291 

l l 

2 3 

0 0 2 ftft  3 9 

. 20. 

32 

17.0  0.0 

. 0. 

i 

0 

3 

0 

0.  0.  11  0 1.0 

0 

080301 

0 4 

1 0 40  74  . 

3135 

4.71 

2 40.  0.  0 

. 0. 

0 

0 

3 

0 

0.  0.  11  0 0.0 

0 

0803  1 1 

J 3 

2 

2 -4074. 3135  14.75 

72.  0.  0 

. 0. 

0 

0 

3 

0 

0.  0.  11  0 0.0 

0 

08032 1 

0 1 

3 

0 0. 

1. 

0 

30.  0.  0 

. 0. 

0 

0 

3 

0 

0.  0.  11  0 0.0 

0 

080331 

0 10  4 

0 0. 

1 . 

00 

25.  1.0 

. 0 . 

0 

0 

0 

0 

o.  o.  u o o.o 

0 

082000 

l 

•90XXY 

PUMP  DESCRIPTION  DATA  CARDS 

0900  1 1 

3 

4 

0 l 0 1868. 

1.0 

4494  3 . 7 

10.  3 

2 2 

200  . 

20000. 

4 7.17 

0. 

0 

09002  1 

3 

4 

0 1 0 1668. 

l .0 

44943.  71Q.2 

•9  6 
4.  - 

200. 

20000. 

4 7.27 

Cl. 

0 

091001 

- 1 1 

0. 

0.  . 1 

0. 

. 

15  .05  .2 

- 

8 

. 3 

.96  .4 

. 98 

. 6 

. 9 7 

091002 

.3 

.9  .9 

. 8 

• 

96  .5  1. 

0 0 

. 

09200 1 

_ "> 

4. 

0. 

0.  1 . 0. 

095011 

0. 

0 

. 0. 

09502  1 

0. 

0 

. 0. 

1 00000 

0 

0 

16  0 

103011 

1 

1 

5 0.0 

1.31 

0. 

25  1.25 

0 .*5 

1 

n 

0 

0.7  5 

1.12 

1.0 

1 . 0 

1 0302  1 

l 

2 

5 0.0 

-0.9 

m 2 

5 -.5 

. 5 

- 

.05 

. 75 

. 4 

1 . 0 

1 .0 

103031 

l 

3 

5-1.0 

2 . 1 

- . 

75  1 . 77  5 

-0  . 

5 

1 

. 54 

-.25 

l . 4 

0. 

1.31 

103041 

1 

4 

5 -1.0 

2 . 1 

- . 

75  1.50 

-0. 

5 

l 

. 1 3 

-.25 

. 92 

0. 

.8  2 5 

103051 

l 

5 

5 0. 

• 4 

n 

5 .53 

m 5 

. 

62 

5 

. 75 

. 80 

1.0 

1 .0 

103061 

1 

s 

5 0. 

. «:  5 

* 

5 . 79 

. 5 

. 

91 

5 

. 75 

. 9 C 

1 .0 

1 . 0 

103071 

l 

7 

5 -1.0 

-i  . 9 

75  -1  . 

-.  5 

- 

. 4 

-.25 

. 1 

0. 

. 4 

103081 

1 

3 

5 -1.0 

-1.8 

. . 

75  -1.72 

-.5 

- 

1 . 

54 

-.25 

-1.2  7 

0. 

- . 9 

103091 

l 

5 0. 

.55 

*» 

5 .63 

. 5 

71 

. 75 

. 86 

1.0 

1 . 0 

103101 

2 

5 0. 

-.55 

■y 

-.  1 75 

• 4 

0 

. 75 

.58 

1.0 

1 . 0 

103111 

3 

5 -I  .0 

1 . 77 

75  1.23 

-.  5 

. 

8 3 

-.25 

. 62 

0. 

.55 

103121 

2 

4 

5 -1.0 

l . 77 

- . 

75  1.52 

-.  5 

1 

. 3 

2 

-.25 

1.13 

0. 

l . 0 

103131 

5 

5 0. 

-.  8 

0. 

2 -.375 

0.  5 

- 

. 1 

. 75 

. 1 5 

1 . 0 

.43 

103141 

■5 

ft 

5 0. 

1.0 

5 . 88 

. 5 

75 

. 75 

. 6 1 

1.0 

.43 

103151 

7 

5 -1  . 

-3  .45 

- . 

75  -2.75 

-.  50 

- 

z . 

0 

— * 2 ^ 

-1.375 

0 . 

- . 8 

133181 

> 

3 

5 -I  . 

-3  .45 

- . 

75  -2.62 

-.  50 

- 

1 . 

58 : 

! -.25 

-1.15 

0. 

55 

• 1 1XXX0 

valve 

DATA  CARDS 

110010 

4 

0 . 

0 0.  0. 

0. 

0. 

110020 

5 

2 

0 0 . 0 

. 0. 

o. 

110030 

-0 

0 . 

0 0 . 0 

. 0. 

o. 

110040 

- 3 

0 

. 0 0.  0. 

0.  0 

• 

• 12XXYY 

LEAK 

TABLE  OA 

TA  CARDS 

120101 

..  i 

9 

30.  0 

1.  1000. 

1 . 

120201 

ft 

5 1 

4.70.  1 . 

2.25 

1 . 

2.5  .667  2. 

75  . 

33 

3 

3 . 

0. 

10000. 

0. 

* 1 3 X XY v FILL  TABLE  DATA  HARDS 


•3XX00 

TRI®  TYPE 

POINTS  ICALC  UNITS 

PORT 

HOP.X 

1 30100 

2 l 

4 1 

LBS/SEC 

1050. 

40  l . 

11324 

1 30101 

0.  364.6676  1.  864. 

66  76  4.0  0 . 

10000 

. 0. 

l 302  00 

1 -6 

4 3 

LBS/SEC 

54  8 . 

. 999 

1 30201 

0.  0.  940 

. 0.  1025. 

-2  76.22  4 6 

3000  . 

-2  76.22 

46 

l 30  300 

7 : 

1 1 3 

CAL  /**IN 

2 12  . 

0.0 

1 30301 

0.  15150. 

59.  14200. 

13  7 . 12  500 

. 17*. 

1 1 550  . 

2 3 7 . 0 4 7 0 . 

1 30302 

268.  7576. 

280.  5682 

. 2 8 7.  3 788 

. 2 9 7. 

1 8 9 4 . 

304.  0.  10000 

l 30400 

7 > 

9 3 

CAL /MI H 

212  . 

0 

. 

1 3040  1 

1 7.6  4 3 55  . 

2 3 . 7 4 30« 

. 26  . 1 42  DI 

. 41.2 

4175. 

42.4  4lft6. 

1 30402 

°8  . 3 708  . 

102  . 36  72  . 

310.  0.  3000.  0. 

• 140C0 

riNTTICS  CONSTANTS 

*14000 

PUR  MUL 

RIL  R EC 

UDU? 

140000  3 l 116.  0.  1. 

•1.1001  SCRAM  TABLE  DATA  CARDS 

141001  -8  4 0.  0.  0.9  0.  1.075  -.4  1.8  -4.  3.7  -SO.  5.7  -23.1  S.  -30.  10.  -30. 

•142001  DENSITY  REACTIVITY  DATA  CARDS 

[42001  -10  . 2 - 1 1 . 5 5 . 3 -8.830  .4  -8.890  . 5 -5.01  . 8 -3.84  .7  -2.51 
l 42002  .8  -1.55  .9  -0.73  l.  0.  1.1  0. 

•143001  DOPPLER  TABLE  DATA  CARDS 

141001  -8  0.  1.85  1000.  0.  2000.  -1.71  3000.  -3.  4000.  -4.28  5000.  -5.43 


t 


TABLE  5 


(cont  . ) 


•UOXXO  R £>\CT  I V l T Y COCFFICISNT  DATA  CARDS 


140010 

148289 

. L4828« 

0. 

0. 

140020 

3* 1 782 

. 341 782 

0 . 

0. 

140030 

019858 

.010858 

0. 

0. 

140040 

34 1 782 

. 34 l 782 

0. 

0. 

1 40050 

148289 

. 1 482  89 

0. 

0. 

• l 5 X X X 

H F AT  SLAB 

DATA 

15001  1 

0 

3 2 0 

0 1 l l 1 

0. 

88  7 . 

42  2 . 

0. 

0. 

0. 

0. 

0. 

0 

l 500  12 

0. 

, 0. 

l 5002  1 

0 

2 2 0 

Oil  ll 

0. 

1 368. 

698  . 

0. 

0. 

0. 

0. 

0. 

0 

l 50022 

0. 

, 0. 

l 5003  l 

0 

4 2 0 0 

l l l l 0. 

1918. 

97  7 . 

0. 

0. 

0. 

0. 

0. 

0 

l 50032 

0. 

. 0. 

l 5004  1 

0 

112  0 

0 l 1 1 l 0 

. 

68  7 . 

385  . 

0. 

0. 

0. 

0. 

0. 

0 

150042  0.  0. 


150051 

41830011 

l l 

1355.  1355.  23° 

. 0. 

0. 

0.  0. 

0 

0. 

150052 

0.  0. 

1 5006  1 

1240011 

1 1 

2298.  2288.  286 

. 0. 

0. 

0.  0. 

0 

0. 

1 50062 

0.  10.42 

150071 

11235001111197.197.  65. 

0* 

0. 

0.  0. 

0 

0. 

1 500  72 

0.  0. 

1 5003 1 

16150011 

1 1 

94  . 94  . 3 1 . 

0. 

0. 

0.  0. 

0 

0. 

150082 

(3.  0. 

150091 

11176001 

1 11  6884.  6884.  95 

. 0. 

0. 

0.  0. 

0 

0. 

150092 

0.  17.1 

150101 

01170011 

1 l 

0 . 7 36  . 46  . 

0. 

0. 

0.  0. 

0 

0. 

150102 

0.  0. 

150111 

2380011 

l 1 

4380.  4380.  184. 

0. 

0. 

0.  0. 

0 

0. 

150112 

0.  0. 

l 50  1 2 l 

0590011 

1 1 

0 450.  36. 

0. 

0. 

0.  0. 

0 

0. 

150122 

0.  0. 

150131 

0993011 

1 1 

0.  450.  36. 

0. 

0. 

0.  0. 

0 

0. 

150132 

0.  0. 

150141 

0 7 9 0 0 1 1 

l 1 

0.  338.  67. 

o. 

0. 

0.  0. 

0 

0. 

150142 

0 . 0. 

150151 

01090011 

I t 

0.  8 38.  67. 

o . 

0. 

0.  0. 

0 

0. 

150152 

o.  o. 

150161 

12  18  1000 

1 l 

11  1096. S 1096.9 

2 5. 

b : 

0.  0. 

0. 

0.  1.375 

1 . 3 

75 

150162 

0.  2.74 

1 501  71 

13181010 

1 1 

1 1 1096.  9 1 096 . 3 

2 5. 

6 7 

0.  0. 

0. 

0.  4.12  5 

4.1 

2 3 

150  172 

2 . 76  5 . 49 

150131 

14  19  10  1 0 

1 l 

11  398.84  398.34 

9.33- 

0 . 0. 

0. 

0 . 6.0  4 

. 0 

150182 

5.51  6.49 

150191 

15131010 

l 1 

l l 1096 . 8 1 096 . 3 

2 5 . 

67 

0 . 0. 

0. 

0.  7.975 

7.9 

75 

150192 

6.51  9.24 

150201 

18191010 

1 1 

11  1096.8  1096.9 

25. 

6 7 

0.  0. 

0. 

0.  10.62 

5 10 

. 62  5 

150202 

9.26  11.99 

150211 

01210011 

1 1 

0.  15175.52  17 

8.03°  0 

. .048 

0 

. . 05  8 8 1 

. 375 

1.375 

1 502  1 2 

0.  0. 

1 502  2 1 

0 1 3 1 1 0 1 1 

l l 

0.  15178.52  17 

8 . 039  0 

. .048 

0 

. .0588  4 

. 12  5 

4.125 

l 50222 

0.  0. 

1 502  3 1 

0 1 4 1 1 0 l I 

1 1 

0 . 55  1 9 . 3 72  64  . 

741 

0 . 

048  0 . 

.0588  6.0 

6.0 

1 502  32 

0.  0. 

150241 

01511011 

l 1 

0.  15175.52  178 

.039 

0. 

. 048 

0. 

.0588  7. 

975 

7.975 

1 50242 

0.  0. 

150251 

01611011 

l 1 

0.  15175.52  178 

.039 

0. 

.043 

0. 

.0588  10 

.62  5 

10.62  5 

1 502  52  0.  0. 


*2  34  56  7890  1 

2 

345678901234567990123456 

7990  1 

2 3*. 

54 

7990  l 2 

34 

567890 

• 1 6000 

C OR  X 

SECTION 

DATA 

♦ 16000 

S LB 

MODE 

CLAD 

OF*  AC 

l 600  l 1 

2 l 

* 

9 15 

0. 

. 1 7 72  1 

160020 

2 2 

- 

9 I 5 

0. 

.26902 

160030 

2 3 

y 

9 1 5 

0. 

.10754 

1 60040 

2 4 

7 

9 1 5 

0. 

.26902 

160050 

25 

7 

9 15 

0. 

.17721 

1 600  1 5 

0 

0 

96.07 

. 3 

.4  .965 

.95 

0 

0 

.05 

0 

0 

0 

.08333 

.9714 

.01 

0 

0 

1 6002  5 

0 

0 

0 . 

. 4 

.6  .96  5 

. 95 

0 

0 

.05 

0 

0 

0 

.08333 

0. 

.01989 

0 

0 

160035 

0 

0 

0. 

. s 

.6  .96  5 

.95 

0 

9 

.05 

0 

0 

0 

.08333 

0. 

. 0 1 9 ? A 

0 

0 

1 60045 

0 

0 

0 . 

.3 

.4  .965 

. 95 

0 

0 

.05 

0 

0 

0 

.08333 

0. 

.01988 

0 

0 

150055 

0 

0 

.1:59 

. 8 

.4  .965 

.95 

0 

0 

.05 

0 

0 

0 

.08333 

.00  12  7 

.01993 

0 

0 

-A' 


— — 


TABLE  5.  (cont.) 


7 7 


* l 7XXYY  SIAR  GEOMETRY  DATA  CARPS 


1 70101 

2 3 

i 

6 0.0  .020071  1 .0 

1 70102 

1 

2 

2 .000721  0. 

1 70 103 

0 

3 

8 .00267  0. 

1 70201 

1 2 

4 

30.0  .4  0. 

1 70202 

0 

5 

1 .0208  0. 

1 70301 

1 l 

5 

10.  .17  0. 

l 7040 l 

1 1 

5 

1 0.  . 12  5 0. 

1 70501 

1 1 

5 

1 0.  .33  0. 

1 70601 

1 1 

5 

10.  .0138  0. 

1 70  701 

1 1 

5 

1 0.  .0625  0. 

1 70801 

1 1 

5 

1 0.  .042  0. 

1 7080! 

1 1 

5 

1 0.  .08  0. 

1 71001 

1 1 

5 

1 0.  .0234  0. 

130100 

-20  500. 

3.341  650.  2.971 

800. 

2.677  950. 

2.439  1100.  2.242 

180101 

1 

2 50 

. 2.078  1400.  1.940 

1 550  . 

1 . 82  3 1 700 

. 1 . 72  4 1 8 50  . 1.639 

1 30102 

■» 

000 

. 1.568  2150.  1.507 

2300. 

1.457  2450 

. 1.415  2600.  1.382 

180103 

3100 

. 1.323  3600.  1.333 

4 100. 

1.406  4600 

. 1.538  51 

00.  1.730 

1 802  OJ 

2 

32.  .41562 

5400  . 

.41562 

1 30300 

1 8 

32.  7.812 

2 12. 

7.992 

392  . 

8.208 

180301 

5 72  . 3.  784 

752  . 

9.540 

9 32  . 

10.404 

1 80302 

1112.  11 .268 

1 2 92  . 

12.492 

1 4 72  . 

13.176 

180303 

1652 . 13.968 

1 3 32  . 

14.  7<ib 

2012  . 

16.128 

1 30304 

2 1 92  . 1 7 . 784 

2 3 72  . 

19.656 

2 5 5 2 . 

21.780 

180305 

2 732  . 24.048 

3092  . 

2 8.  908 

3 3 60  . 

33.120 

1 80401 

-5  32 

. 30.  212.  29.5  392.  28.3  572.  26. 

6 752.  24.7 

130501 

zoo. 

3.  33  1200  . 12.92 

1 90100 

1 6 

32  . 34.45 

1 

22.  38.35 

212. 

40.95 

190101 

392.  43.55 

7 

52 . 46.8 

2012  . 

51.35 

190102 

2 7 32  . 52.85 

3092  . 

56.55  3452 

. 83.05 

190103 

38  12  . 72  .8 

4352.  $9.7 

4 5 32  . 

94.25 

180104 

4712.  98.15 

*.892  . 100.1 

5 144. 

10  1.4 

190105 

8000.  101.4 

1 902  00 

2 

32  . .0000  75 

5400  . 

. 0000  75 

1 90300 

5 

0.0 

2 8 . 392 

190301 

1 480 

. 3 

34.478  1875.0  85.176 

190302 

1 '9  7 

. 5 

34.476  3500.0  34.476 

190401 

-7  130 

. 56.9  350.  80.3  450.  52.3 

5 30  . 6 5 

.2  620.  67. 

2 710.  70.: 

800.  77.5 

19050! 

-10  68 

, 5 

: . 9 200  . 56  . 7 400  . 6 1 .6 

600  . 64 

. 800.  66. 

1000.  .87 

1 9o  5 0 2 

12 

00. 

88.4  1400.  71.3  1800.  75.8  1800.  90.6 

200101 

-2  0. 

3 • 71 8E-6  5000 . 1 . 2 6 5 3E-5 

200201 

-2  0. 

0. 

5000.  0. 

200301 

-4  0 

. 

3.094E-6  1852.  4.706E- 

6 1653 

. 5.389E-6 

5000  . 5 . 38  0 E-6 

200401 

-2  0. 

0 

. 5000.  0. 

200501 

-2  C . 

0 

. 5000.  0. 

250001 

1 2 3 

3 . 

06  1 5 

260001 

1 2 3 

4 6 

666111E7110000 

0 0 0 

0 0 

2 7000  l 

- 1 - : 500. 

0.  22  00  . 98.7  1820.  197 

. 1 730  . 

395.  1660. 

592.  1600. 

*<30  . 

2 70002 

1 5 4 0. 

^8  7 

. 1480.  1134.  1440.  1382 

. 1400  . 

1579.  1370 

. 1777.  1335 

. 1**4. 

2 70003 

13  10. 

: i 

72  . 12  90  . 2 3 6 9 . 

2 3000  l 

-4  3 4 

. 

0.  34.  150.  21.  200. 

2 l . 

1000  . 

2 90001 

-14  2500. 

0.  2200.  99.7  1820.  197 

. 1 730  . 

395.  I860. 

592.  1600. 

790  . 

2 90  0 02 

1 540. 

98  7 

. 1480.  1184.  14.0.  1382 

. 1400. 

1579.  1370 

. 1 7 7 7 . 1 3 3 5 

. 1*74. 

2 90003 

1310. 

2 l 

72  . 12  90  . 2 369. 

30000 1 

-14  131. 

0.  131.  98.7  62.  197.  47 

. 3 9 5 . 

89  . 592  . 0 5 

. 790 . 116. 

98  7 . 

300002 

132. 

1 1 94  . 1 3 7.  1 382  . 1 3 1 . 1 5 79 

. 12  1. 

1 777.  100. 

1924.  76.  2! 

72  . 

300003 

59.  2 

369 

• 

TABLE  6.  STANDARD  MODEL  INPUT  DATA  DECK 


7« 


-BUR  BEST  ESTIMATE  2005RRK  CRITICAL  FLOW  - SONIC 
*2  3456  789012  34  56  78  90  II  3456  789012  34  56  7890  12  3456  7S90 

*10001  LDM  ED  I NTS  TRP  VOL  BUB  TD  V JUN  PMP  CKV  NLK  FLL  SLBV  COM  MAT  COP. 
010001  -2  9 5 10  23  4 0 3 3 2 4 2 4 2 5 1 0 5 5 0 0 

010002  3388.  1.0 

020000  AP  11  JW  27  JW  28  JW  29  JW  16  ML  4 SR  23  CR  23  TD  27 

030010  1 100  5 0 .001  .00001  .1  2000. 

030020  5330  .01  .000001  l. 

030030  20  10  2 0 .01  .000001  80. 

030040  20  25  1 0 .01  .000001  130. 

030050  20  25  l 0 .01  .00001  2000. 

030003  50 

040010  1100  60.0  0. 


040020 

l 

-4 

3 0 40.0  0 

040030 

2 

1 

0 0 0.0  0 

040040 

3 

1 

0 0 0.001  0 

040050 

4 

1 

0 0 .002  0 

040060 

5 

1 

0 0 1.0  0 

0400  70 

6 

1 

D 0 1 . E6  0 

0 4 0 0?  0 

7 

-5 

4 0 21.5  26 

. 5 

040090 

3 

l 

0 0 .001  1 

. E6 

040 1 00 

9 

-5 

4 0 21.3?  1 

2 0 . 

• 5 r x x x 

VOLL' 

ME  DATA  CARDS 

0500  1 l 

0 

1039.  -1 . .091632 

1 54  8 . 

22.25 

22.25 

0 

^ >6  ' 

0. 

29.92 

0 

■■  521 

2 

0 

1023.23  -1.  .002448 

6 703  . 

12.25 

5.2  7 

0 

497. 

0. 

41  . 75 

0 

050031 

1 

0 

102  4.98  - 1 . . 999 

3 7 7 7 . 

21.25 

2 1.25 

0 

497. 

0. 

51.25 

0 

050041 

3 

0 

1 03  1.93  532.3  -1  . 

3177. 

32.23 

32.23 

0 

165. 

0. 

10.02 

0 

050051 

0 

0 

102  4 . 5 32  . 3 -l . 

125.25 

26. 

26. 

0 

5.24 

0. 

- 1 4 . 

0 

05006  L 

3 

0 

1126.17  533.25  -1 . 

83. 

3. 

3. 

0 

4 . 0 

0 . 

-14. 

0 

050071 

0 

0 

1226.17  533.25  -1. 

366  . 

40.53 

40.53 

0 

3.2 

0. 

-14. 

0 

0 5 5081 

0 

0 

1024 . 532.30  -1 . 

136. 

2 6 . 

26. 

0 

3.67 

0. 

-14. 

0 

050091 

0 

0 

1126.17  533.25  -l. 

83  . 

3. 

3 B 

0 

4 . 0 

0 . 

-14. 

0 

C 5 0 1 0 1 

0 

0 

1226.17  533.25  -l. 

366  . 

40.53 

40.53 

0 

3.2 

0 . 

-14. 

0 

050111 

■> 

0 

1062.18  532.53  -1 . 

2 131.5 

17.20 

17.20 

0 

120. 

0 . 

0. 

23 

u 5 0 1 2 l 

2 

0 

1054. 87  549.337  -1. 

22  3 . 

2.85  2 

.35  0 

31 

.091 

. 04  73 

1 7.3  2 

1 3 

050  1 3 l 

0 

1 052.9  7 - 1 . . C3  7 1 5 7 

2 2 3 . 

2.75  2 

. 75  0 

8 l 

. 09  1 

. 04  73 

20.67 

1 4 

050141 

0 

1051.52  -1.  .050567 

9 1.09 

1 . 1 

0 

31 

. 09  l 

. 04  73 

2 3.42 

l 5 

050  1 5 1 

0 

1049.97  -1.  .091994 

2 2 3 . 

2.752 

. 75  0 

31 

. 09  l 

. 04  73 

2 4.42 

l 6 

050161 

0 

1 04  7.85  - 1 . . 139875 

2 2 3 . 

2.752 

. 75  0 

31 

. 09  1 

. 04  73 

2 7 « 7 

1 

350171 

» 

0 

1 05  . . 2 8 5 32  . 5 3 - 1 . 

442.47 

18.12 

18.12 

C 

52.30  1 0. 

0. 

1 9 

050131 

2 

0 

1045.96  532.87  -1. 

534.53 

12.35 

12.35 

0 

43.63 

5 0 . 

17.92 

1 

050191 

0 

0 

1059.2  532.53  -1 . 

115. 

15.7 

15.7 

0 

1 0 . 69 

0. 

10.02 

0 

05020 1 

0 

3 

1059.2  532.53  -l. 

115. 

15.7  15.7 

0 

1 9 . 69 

0. 

10.02 

0 

0502  l l 

n 

0 

1030.9  532.30  -1 . 

10.  75 

1 1 

6*4 

■>  # 2 

0 

5.24 

0. 

10. 

0 

0 5 0 2 2 1 

4 

o 

14.7  100.  0.6 

3 . 4 3 E 6 

2 50  . 

0.0 

1 

1000. 

100. 

-10. 

0 

0 502  3 l 

0 

1055.48  532.53  -1 . 

20  . 5 

1.25 

l . : 5 

0 

3 1.09 

0. 

16.77 

12 

•6000  LIOUID  LEVEL  CALCULATION  CARD 
C60C00  11  23  12  13  14  15  18  17  [8  19  20  l 
380001  0.0  0.0  0.0  0.0  0.0 

*80002  WALLIS  SLIP  CORRELATION  CARD 
•80002  UALSC1  VALSC2 
*8  XXX 1 BUBBLE  DATA  CARDS 
•8XXX1  ALP'!  VBUS 
06001  1 1.0  0.0 
060021  0.0  3.0 

060031  0.3  3.0 

0600-*  l 0.  1.E6 

*60001  SLIP  VELOCITY  CARD 
•8XXXY  JUNCTION  DATA  CARD 

•23456  78  901234567?  90  123456789012345679  90  123456789012345678  90  1234567S9012 


0 3 00  1 l 

l 

0 

0 

29583. 

191  . 

1 52.1 

7 9.  0.  0. 

1 

-1 

3 

0 

0.0 

0.  1 

0 l . 0 0 

C8002  1 

y 

3 

0 

0 

4074.3135  312.  52 

.25  >1.  0. 

0 . 

l 

-l 

3 

0 0 

0 0. 

l 0 0.0 

030031 

4 

0 

0 

25444. 

165  . 

42 .0 

0 . 0 . 0 . 

1 

-1 

3 

0 

0.0 

0.  1 

0 0.0  0 

0300.4  i 

4 

20 

0 

0 

10041.5 

3.11 

2 2 5 . 

30.  .2084 

1 . 

l 7 

2 

-l 

0 2 

0.0 

100. 

0 9 0 0 5 1 

4 

2 1 

0 

0 

4 750* 

3.8  7 

11.0 

0.  0.  0. 

1 

-1 

3 

0 

0.0 

0.  1 

0 0.0  0 

C 8 00 6 l 

5 

*) 

-l 

0 

4 750  . 

3.87 

-13. 

0.  0.  0. 

0 

-I 

3 

0 

0.0 

0.  1 

0 0.0  0 

: 3oo 7 1 

6 

7 

l 

0 

4 750  . 

3.2 

-1  3. 

0.  0.  0. 

0 

-1 

3 

0 

0.0 

0.  1 

0 0.0  0 

7 

20 

0 

0 

4 750  . 

.539 

2 5.  4 

7 . . 2 3 73  6 

. a 

0 

-l 

0 

: 0 

0. 

0 0.0  0 

0 8006  t 

4 

9 

0 

0 

4 7 5 0 . 

3.67 

11.0 

0.  0.  0. 

l 

-1 

3 

0 

0.0 

0.  1 

0 0.0  0 

080!  H 

8 

9 

_ 1 

0 

4 750  . 

3.6  7 

-13. 

0.  0.  0. 

0 

- 

3 

0 

0. 

0.  1 

0 0.00 

0801 1 1 

9 

10 

2 

0 

4 750  . 

3.2 

-13. 

0.  0.  3. 

0 

-1 

3 

0 

c . 

0.  i 

0 0.00 

080121 

10 

: 9 

0 

0 

4750. 

.538 

2 5 . 4 7 . . 2 3 73  6 

. PI 

3 

-1 

0 

2 0 

A , 

0 O.0  0 

080131 

4 

19 

0 

0 

10041.5 

3.11 

2 2 5 . 

30.  .2084 

t . 

1 7 

2 

- 1 

0 : 

0.  0 

1 0 0. 

080141 

: 3 

12 

0 

0 

26639. 

8 1 . 

09  1 7 

.92  0.  0. 

0. 

n 

-1 

) 

0 0 

0 • 

0 1.0  0 

080151 

• ^ 

! 3 

0 

0 

26639. 

81  . 

09  20 

.87  0.  0. 

0. 

0 

-1 

3 

0 0 

0 . 

101.00 

080161 

1 3 

1 4 

0 

0 

2 66  3°. 

3 1 . 

09  2 3 

.42  0.  0. 

0. 

0 

-1 

3 

0 0 

ft 
'j  • 

0 1 .0  0 

i 


' 

7 o 

TABLE  6.  (cont.) 


080 1 71 

1 4 

15 

0 

0 

2*6  39 

. 81.09 

2-. 42 

0 . 

0. 

0. 

0 

- 1 

3 

0 

0. 

0. 

1 

0 

1 . 

. 0 

0 

080131 

15 

1 6 

0 

0 

2 66  39 

. 81.09 

2 7.17 

0 . 

0. 

0 

-I 

3 

0 

0. 

0. 

I 

0 

1 . 

.0 

0 

080191 

16 

1 

0 

0 

2 6639 

. 37.13 

2 9.92 

0. 

0. 

0 . 

0 

-1 

3 

0 

0 . 

0. 

i 

0 

1 . 

. 0 

0 

080201 

1 1 

l 7 

0 

0 

2 944  . 

2.246 

17.0 

0. 

0. 

0. 

1 

-1 

3 

0 

L . 

0 0. 

1 0 ( 

).0 

080211 

1 7 

18 

0 

0 

2 944  . 

26.20 

18.02 

0. 

0. 

0. 

0 

-1 

3 

0 

0. 

0. 

l 

0 

0. 

. 0 

0 

080221 

18 

1 

0 

0 

2 944  . 

*3.635 

30.02 

0. 

0. 

0. 

0 

-l 

3 

0 

0 . 

0. 

1 

0 

1 . 

. 0 

0 

080231 

19 

11 

0 

0 

1 4 79  1 

.5  19.69 

10.12 

0. 

0. 

0. 

l 

-l 

3 

0 

0 . 

0. 

1 

0 

0 , 

. 0 

0 

0802-1 

20 

1 1 

0 

0 

l 4 79  l 

.5  19.69 

10.12 

0. 

0. 

0 . 

1 

-l 

3 

0 

0 . 

0. 

1 

0 

0 , 

. 0 

0 

080251 

2 1 

5 

0 

1 

4 750  . 

3.67 

11.0 

0. 

0. 

0. 

0 

-1 

3 

0 

0 . 

0. 

l 

0 

0. 

.0 

0 

080261 

3 2 

i 

0 

3 

0. 

.484 

55.0 

1 . 

1 . 

1 . 

0 

-1 

0 

0 

0. 

0. 

1 

0 

0. 

. 0 

0 

0802  71 

21 

22 

0 

4 

0. 

3.67 

1 1 . 

0. 

1 . 

. 5 

0 

2 

2 0 0 

1 . 

• 9 _i_ 

0 0.< 

) 0 

080281 

22 

5 

0 

4 

0 . 

3.67 

1 1 . 

0. 

. 5 

l . 

0 ‘ 

T 

2 0 0 . 

* 9 1 

0 0.( 

3 0 

080291 

1 1 

23 

0 

0 

2 66  39 

. 20.32 

1 7.0 

0 . 

0. 

0. 

r 

3 

0 

0. 

0. 

l 

0 

1 . 

. 0 

0 

080301 

0 4 

1 

0 

40  74  . 

3135  4.712 

*0. 

0. 

0. 

0. 

0 

-1 

3 

0 

0. 

0. 

l 

0 

0 , 

.0 

0 

08031 1 

0 3 

2 

2 

- 

40  74.  3 1 3 5 1 4 . 75 

72  . 

0. 

0. 

0. 

0 

-1 

3 

0 

0 . 

0. 

l 

0 

0 

.0 

0 

080321 

0 l 

3 

0 

0 

. 

1.0 

30. 

0. 

0. 

0. 

0 

-l 

3 

0 

0. 

0. 

l 

0 

0 , 

.0 

0 

080  331  0 10  4 0 0.  1.00  25  . 1 . 0 . 0 . 0 - 1 0 0 0 . 0 . 1 0 0.0  0 

•90XXY  PUMP  DESCRIPTION  DATA  CARDS 


0900  l 1 

3 

4 0 10 

1 668  . 

1 . 0 4494  3 . 710.3 

22 

200  . 

20000. 

4 7.17 

0. 

0 

09002 1 

3 

4 0 10 

1 668. 

1.0  44943.  710.3 

*>  y 

2 00  . 

20000 . 

47.17 

0. 

0 

061001 

-1  I 

0.  0. 

1 0. 

. 1 5 

. 05  . 2 4 . ’ 

. 3 

.96  .4 

. 98 

. 6 

. 6 7 

091002 

.8  .9 

9 . 8 

. °D 

.5  1.0  0. 

092001 

-2 

0.  0.  1.  0. 

065011 

0 . 

0.  0. 

095021 

0. 

0.  0. 

1 00000 

0 

0 16  0 

10301  1 

1 

l 5 0.0 

1.31 

0.25 

1.25  0.5 

l 

.20 

0. 75 

1.12 

l .0 

l . 0 

1 0302  1 

l 

1 5 3.0 

-0 . 9 

.25 

-.5  .5 

- 

.05 

. 75 

. 4 

I . 0 

1 . 0 

103031 

l 

3 5-1.0 

2 . 1 

-.  75 

1.775  -0.5 

1 

. 54 

-.25 

l . 4 

0. 

1.31 

1030*1 

l 

4 5-1.0 

2 . 1 

-.  75 

1.50  -0.5 

l 

. 1 3 

-.25 

. 92 

0. 

.82  5 

103051 

1 

5 5 0 . 

. 4 

. 2 5 

.53  .5 

. 

*2  5 

. 75 

. SO 

i . 0 

1 . 0 

103U81 

1 

6 5 0. 

. 82  5 

.25 

.79  .5 

. 

815 

. 75 

.90 

1.0 

1.0 

1 030  7 1 

l 

7 5-1.0 

-1.8 

-.  75 

- 1 . -.5 

- 

. 4 

-.2  5 

. 1 

0. 

. 4 

103081 

1 

8 5-1.0 

-1.8 

-.  75 

- 1 . 72  -.5 

- 

1 . 54 

-.25 

-1.27 

0. 

-.  9 

103091 

2 

1 5 0. 

.55 

.25 

.63  .5 

. 

71 

. 75 

. 86 

1 . 0 

1 . 0 

l U 3 1 0 l 

2 

2 5 0 . 

-.55 

• 2 

-.175  .4 

. 

0 

. 75 

.58 

l .0 

1 .0 

103111 

y 

L 

3 5-1.0 

l . 77 

-.  75 

1.23  -.5 

. 

83 

-.25 

. 62 

0. 

.55 

103121 

y 

4 5-1.0 

1.77 

-.  75 

1.5:  -.5 

I 

. 32 

-.25 

1.13 

0. 

I .0 

103131 

•y 

5 5 0 . 

- . 8 

0.2 

-.  3 75  0 .5 

- 

. 1 

. 75 

. 1 5 

1.0 

.43 

103141 

y 

5 5 0. 

1 . 0 

. 2 5 

.88  .5 

. 

75 

. 75 

. 6 l 

l . 0 

.43 

103151 

y 

& 

7 5 -1  . 

-3.45 

-.  75 

-2.75  -.50 

. 

2 . 0 

-.2  5 

-1.375 

0. 

-.  9 

103141 

y 

8 5 -l  . 

-3.45 

-.  15 

-2.62  -.50 

- 

1.583 

-.25 

-1.15 

3. 

-.55 

* ! ! xx.ro 

valve  data  cards 

11001  0 

0.  0 0.  0 

. 0. 

0. 

110020 

5 

2 0 0 . 

0.  0. 

0. 

1 100  30 

-0 

C.  0 0. 

0.  0. 

0. 

1 100-0 

- 3 

0.  0 0.  0 

. 0.  0. 

• 12XXYY 

LEAP  TABLE  DATA  CARDS 

120101 

-2 

9 30.  C 

i.  1000.  1. 

120201 

6 

5 14.70.  1 

. 2.25 

1 . 2.5 

.867  2.75  .33 

n 

3.  0. 

1 oooc . 

j • 

* i oxx*" 

F 

ILL  7A3LE  DATA  OAR 

ns 

* 3XX00 

7RI®  TYPE  PO  IN 

TS  ICAL 

C I'N  ITS 

PORT  HORX 

1 30 100 

2 

1 4 

l 

LBS/SEC  1050.  4 

C 1 

.1132 

4 

1 30101 

0. 

96  4 .-56  76  1 

. 864.6676  4.0 

0.  10000.  0. 

1 30200 

y 

2 4 

3 

LBS/SEC  548. 

94  9 

130201 

0. 

0.  940.  0. 

102  5 . 

-2  76  . 2 2 

46  3000  . -2  76 

y 

2 4 6 

1 30  300 

7 

2 1 1 

3 

CAL/M 

TN  212.  0. 

0 

1 30301 

0. 

15150.  59. 

14200. 

137.  12500.  178.  11550 

. 237 

. 9*70. 

1 30302 

268 

. 7576.  280 

. 5682  . 

287.  3788.  297.  189*. 

30* . 

0.  10000.  0. 

1 30-00 

* 

2 9 

3 

CAL  /M 

IN  212. 

0. 

1 30401 

17. 

6 4355. 

7 4309. 

2 6.1  4 

291 . 41.2  417 

5 . 

42  . 4 

4 166. 

l 30*02 

68  . 

3709.  102. 

36  72  . 

310.  0. 

3000.  0. 

• 1*000 

Kill 

ETICS  CONST 

ANTS 

•14000 

PUR 

MUL  S I L 

REC 

UDUF 

1 40000 

3 

0 i 1 6. 

0. 

l . 

•141001 

SCRAM  TA3LE  ^ATA  CARDS 

141001 

-8 

* 0.  0.  D . 

9 0.  1 . 

075  -.4 

1.8  -* . 2.7 

- 1 

C . 5 . 

7 -29.1 

8 . -30. 

. 10. 

-30. 

*14200  1 DENS  T T V REACTIVITY  DATA  CARDS 

1 *200  1 -1  0 . 2 -1  1 . 5 5 . 3 -8.830  . 4 -4.880  .5  -5.I'I  .4  -3.44  .7  -2.51 
142002  .8  -1.55  .9  -0.73  l.  0.  1.1  0. 

*141001  DOPPLER  TARIE  DATA  CARDS 

1*3001  -4  0.  1.85  1000.  0.  2000.  -1.71  2000.  -3.  *000.  -4.28  5000.  -5.43 


TABLE  6 


(cont.  ) 


M4QXX0  REACTIVITY  COEFFICIENT  DATA  CARDS 


I 400  l 0 

.148289  .148289 

0.  0. 

140020 

. 34  1 782  . 34  1 782 

0.  0. 

140030 

.019858  .019858 

0.  0. 

140040 

. 34  1 782  . 34  1 782 

0.  0. 

140050 

.148289  .148289 

0.  0. 

*1 5XXX 

HEAT  SLAB  DATA 

15001  1 

0 3 200000 

. 63  7 . 

42  2 . 

0. 

0 . 

0. 

0. 

0. 

1 500  1 2 

0.  0. 

15002  1 

0 2 200000 

. 1368. 

698  . 

0. 

0. 

0. 

0. 

0. 

1 50022 

0.  0. 

150031 

04200000. 

1916. 

97  7 . 

0. 

0. 

0. 

0. 

0. 

150032 

0.  0. 

15004  1 

0 11  2 00000. 

68  7 . 

385  . 

0. 

0. 

0. 

0. 

0. 

1 50042 

0.  0. 

150051 

4 13  30000  1355.  1355. 

2 39 

. 0. 

0. 

0. 

0. 

0. 

150052 

0.  0. 

1 5006  1 

1 2 4 J 0 0 0 22S8.  2288. 

286 

. 0. 

0. 

0. 

0. 

0. 

1 50062 

0.  10.42 

1 500  71 

11  23  5 0 0 0 0 197.  197. 

65. 

0. 

0. 

0. 

0. 

0. 

1 500  72 

0.  0. 

l 50031 

16  l 5000094 

. 94 . 

3 l . 

0. 

0. 

0. 

0. 

0. 

1 50082 

0.  0. 

1 50091 

U 1 7 0 00006884.  6884 

. 95 

. 0. 

0. 

0. 

0. 

0. 

1 50092 

0.  17.1 

150101 

0 11  700000. 

736.  46. 

0. 

0. 

0. 

0. 

0. 

150102 

0.  0. 

150111 

2 380000  4380.  4380. 

184  . 

0. 

0. 

0. 

. 0. 

0. 

150112 

0.  0. 

150121 

05900000 

450  . 

36  . 

0. 

0. 

0. 

0. 

0. 

150122 

0.  0. 

150131 

08900000. 

45U. 

36. 

0. 

0. 

0. 

0. 

0. 

1 50132 

0 . 0 . 

150141 

0 7900000. 

83  8 . 

6 7. 

0. 

0. 

0. 

0. 

0. 

150142 

0.  0. 

150151 

0 10  9 0 0 0 0 0. 

838  . 

67. 

0. 

0. 

0 . 

0. 

0. 

150152 

0.  0. 

150161 

1218100000 

1096.8  1096.3 

2 5 . 

6 7 

0 . 

0. 

0. 

150162 

0.  2.7, 

1501?! 

13  l .3  10  1 0 0 0 

1096.3  1096.3 

2 5 . 

6 7 

0 . 

3 . 

0. 

150172 

2 . 76  5.4  9 

150181 

14  13  10  1 000 

398.84  398.94 

9.334 

0. 

0. 

0. 

150182 

5.51  6.49 

150191 

15  13  10  1 0 0 0 

1096.8  1096.8 

25. 

67 

0. 

0. 

0. 

150  172 

6.51  9.24 

150201 

16  IS  10  l 0 0 0 

1096.3  1096.8 

25. 

67 

0. 

0. 

0. 

1 50202 

9.26  l 1 . 99 

150211 

0 12  1 0 0 0 0 0. 

15175.52 

178.039  0 

. 

048 

0 . 

150212 

0.  0. 

1502  2 1 

0 l 3 1 1 0 0 0 0. 

15175.52 

1 78.039  0 

. , 

043 

0. 

150222 

0.  0. 

150231 

0 14  1 1 0 0 0 0 . 

55  1 3 . 3 72 

64  . 

741 

0 . 

048 

0. 

.0 

1 502  32 

0.  0. 

150241 

0 15  1 1 0 0 0 0. 

15  175.52 

I 78 

.030 

0. 

.043 

0 . 

150242 

0.  0. 

150251 

0 16  1 10000. 

15175.52 

l 78 

.039 

0. 

.048 

0 . 

150252 

0.  0. 

*2345673  90  12345673  90  12345678901234567 

9901 

2 34  56  7390  12  3 

0. 
0. 
0. 
0 . 
n . 
0. 
0. 
0. 
0. 
0. 
0. 


0.  4.12  5 
0.  5.0  6 . 
0.  7.375 

0.  10.52! 

.0538  1, 

.0538  4« 
588  6.0  ( 
.0588  7.1 


• 1 5000  CORE  SECTION  DATA 


• 1 6000 

S L 3 

NODE 

CLAD 

OFRAC 

16001  0 

•9  ] 

7 

9 

15 

0. 

. 1 772  1 

160020 

22 

7 

9 

l 5 

0. 

. 2 6 9C  2 

160030 

2 3 

7 

9 

l 5 

0. 

. 1 0 754 

160040 

•>  ^ 

7 

9 

15 

0. 

.26902 

1 600  50 

2 5 

+ 

9 

15 

0. 

. 1 772  1 

1.375 
4.125 
0 

7.375 
10.625 

375  1.375 
125  4.125 

.0 

75  7.375 
625  10.625 


TABLE  6.  ( c on C • ) 


• 1 7XXYY 

SLAB  GEOMETRY  DATA  CARDS 

1 70101 

2 3 1 

6 0.0 

.0200  71  1.0 

170102 

1 2 

T 

.00072  1 0. 

170103 

0 3 

R 

.00267  0. 

170201 

1 2 4 

3 0.0 

.4  0. 

1 702  02 

0 5 

1 

.0208  0. 

1 70301 

1 1 5 

1 0. 

.17  0. 

1 70401 

1 1 5 

1 0. 

.125  0. 

1 70501 

1 1 5 

1 0. 

. 33  0. 

1 7060  1 

l l 5 

1 0. 

.0138  0 . 

1 70  70  1 

I 1 5 

1 0. 

.062  5 0. 

1 7080  1 

l 1 5 

1 0. 

.042  0. 

1 7090  1 

1 1 5 

1 0. 

.08  0. 

171001 

l 1 5 

t 0. 

.0234  0. 

180100 

-20  500. 

3.341 

650.  2.971 

300.  2 . 

6 7 7 9 5 0. 

2.439  1100.  2. 

242 

1 90101 

12  50  . 

2.078 

1400.  1.940 

1550.  1. 

82  3 1 700 

. 1 . 72  4 1 850  . 1 . 

639 

130102 

2000  . 

1 . 568 

2150.  1.507 

2300.  1. 

4 5 7 2 4 50 

1.415  2600.  1. 

382 

180103 

3100. 

1 . 32  3 

3600.  1.333 

4100.  1. 

406  4600 

1.538  5100.  1. 

730 

1 80200 

2 

32  . 

.41562 

5400  . 

. 4 1 562 

1 80300 

IS 

32  . 

7.812 

2 12. 

7.992 

392.  3.208 

180301 

5 72  . 

8 . 784 

75  2 . 

9 . 540 

932.  10.404 

1 80302 

1112 

11.263 

1 2 92  . 

12.492 

1 4 72  . 1 3 . 1 76 

180303 

1652 

. 13.968 

1 8 32  . 

14. 796 

2012  . 16.123 

1 80304 

2 I.  92 

. 17.784 

2 3 72  . 

19.656 

2552 . 21.780 

1 80305 

2 7 32 

. 24.048 

3092  . 

2 8 . 908 

3360.  33.120 

1S0401 

-5  32.  30 

. 212.  2 

9.5  392 . 28.3 

572.  26.8 

752 . 24.7 

l 8050  1 

-2  200.  3 

.33  1200 

. 12.92 

190100 

1 6 

32.  34.45 

122 

. 38.35 

212.  40.95 

180101 

392.  43.55 

752 

. 4 6.9 

2012.  51.35 

190132 

2 732  . 52  . 65 

3092.  56 

.55  3452 

. 63.05 

1 90103 

38  1 2 . 72  . » 

4 3 5 2 

. 8°  . 7 

4532.  94.25 

190134 

4712.  98.15 

4 8 °2 

. 100.1 

5144.  101.4 

190105 

8000 . 101.4 

l 902  00 

2 

32  . 

. 0000  75 

5400  . 

.0000  75 

1 90300 

5 0.0 

2 8 . 392 

19030  1 

1 48O  . 3 

34.4  76 

1675.0  35. 

1 76 

1 90302 

1787.5 

34.476 

3500.0  34. 

4 76 

1 3040  1 

-7  130.  56 

.9  350. 

60.3  450.  62.3 

530.  65.2 

620.  67. 

2 710.  70.2  500.  77 

. 5 

1 9050  1 

-10  63.  52 

. 8 200  . 

56.7  400 . 61.8 

600.  64. 

800.  66. 

1000.  .67 

l 90502 

12  00  . 

68.4  1 400 . 71.8  1600 . 

75.8  1300. 

80 . 6 

2 00101 

-2  0.  3.713E-6  5000.  1.2653E-5 

2 002  0 1 

-2  0.  0.  5000.  0. 

200301 

-4  0.  3 

. 0 9 4 C - 6 

1552.  4.7C6E 

-6  1653. 

5 . 389E-6 

5000.  5.389E-6 

20040  1 

-2  0.  0. 

5000.  0 

. 

300501 

•2  0.  0. 

50C0.  0 

- 

APPENDIX  B 


PLOT/RESTART  DATA  RETRIEVAL  PROCRAM 


This  appendix  describes  a Fortran  IV  program  written  to 
extract  data  from  the  RE  LAP  4 Plot/Restart  tape.  During 
RELAP4  problem  execution  there  are  two  types  of  information 
written  to  the  Plot/Restart  tape.  They  are  called  common 
block  records,  and  plot  records.  Common  block  records 
contain  information  required  to  restart  the  problem.  Plot 
records  contain  data  that  may  be  useful  for  graphical 
display  of  the  transient  parameters  of  interest. 

A plot  record  is  written  to  the  tape  at  time  intervals 
specified  by  the  user  (ie.  each  time  a minor  edit  is 
printed,  a plot  record  is  written  to  the  Plot/Restart  tape). 
The  length  of  a plot  record  depends  of  the  system  being 
modeled  and  can  be  calculated  as  follows, 

LEN  - 21  + 20  2 4 (N  VO  L ) + 16(NJUN)  a-  20(NSLP) 

where. 


LEN  - length 

0 f 

plot  record 

i n 

computer  wo  r d s , 

NVOL  « number 

0 f 

volumes  in 

the 

R EL A P 4 rod  e 1 , 

NJUN  - number 

o f 

Junctions 

in  the  RELAP4  model. 

N S LB  • numb e r 

o f 

heat  slabs 

i n 

the  RFLAP4  model 

P 3 


For  the  problem  used  in  this  study  (ie.  uope  Creek) 
NVOL  ■ 2 3,  N’JUN  « 33,  NSLB  ■ 25,  and  the  the  plot  record 
length  (LEN)  is  thus  1621  computer  words.  A computer  word 
is  four  bytes  or  32  alpha-numeric  characters.  Future 
reference  to  word  will  imply  computer  word,  unless  otherwise 
indicated . 

The  following  is  an  explanation  of  the  reason  for  the 
numbers  used  in  the  above  equation;  The  first  word  of  each 
plot  record  contains  the  letters  'PLOT'.  The  next  twenty 
words  of  the  plot  record  contain  the  problem  title  and  the 
number  of  junctions  and  volumes  in  the  system.  The  next 
twenty  words  are  twenty  system  parameters,  such  as 


normalized 

power,  total  energy 

stored  in 

the 

fuel, 

transient 

time,  etc. 

(These  parameters 

are  listed 

o n 

T T -2  2 8 of 

the  PELAP4  Users  Manual).  The  next  24  X ’.’VOL  words  (in  this 
case  552  words)  contain  the  24  volume  parameters  listed  on 
page  11-229  of  the  Users  Manual.  This  information  is 
written  to  the  tape  as  follows.  Words  1-23  are  the  average 
pressures  in  volumes  1-23,  words  24-46  are  the  total  mass  in 
volumes  1-23,  etc.  The  next  16  X NJl’N  words  (in  this  case 
528  words)  contain  the  16  junction  parameters  listed  on  page 
11-230  of  the  Users  Manual.  This  information  is  written  in 
the  same  order  as  the  volume  information  (le.  words  1-33  are 
the  mass  flow  rates  for  junctions  1-33,  etc).  The  next  20  X 
NSLB  words  (500  words  in  this  case)  contain  the  twenty  heat 
slab  parameters  listed  on  pages  11-230  and  11-231  of  the 


L i --- i 


r 


Users  Manual*  These  parameters  are  written  in  the  same 
order  as  the  volume  and  junction  data. 

Table  7 is  a listing  of  the  program  written  to  extract 
plotting  data  from  the  Plot/Restart  tape.  The  numbers  In 
the  left  column  are  sequence  numbers  used  in  the  folio  v/ing 
program  description.  The  program  utilizes  a Penn  State 
library  subroutine  called  PRECL/FRUF  (21). 

The  name  of  the  Plot/Restart  tape  from  which  data  is  to 
be  extracted  must  be  entered  on  lines  55  and  1250  (in  the 
example  shown  the  tape  name  is  EDIT2).  The  data  set  name 
(DSN)  on  line  1250  must  he  the  same  as  the  as  the  P S N used 
to  write  information  to  the  tape  during  problem  execution 
(in  the  example  DSN «T POUT). 

The  first  step  of  the  program  (line  400)  is  the  CAL!. 
PRECL  statement.  This  statement  loads  the  entire  first 
record  from  the  Plot/Restart  tape  into  a logical  a r r a v 
called  REC.  The  length  of  the  record  is  also  returned 
the  variable  called  L E N . A 00  loop  is  then  entered 
450)  where  the  following  steps  are  taken.  line  ^ * 
the  first  four  characters  of  the  arrav  or  . 
characters  are  not  the  letters  'PLOT'  then  r-e  r 
a plot  record,  and  execution  skips  to  line  - 
CALL  FBUF  statement-  This  statement  1 a 
of  the  tape  into  the  array  REC  and 
450.  This  process  is  continued  unt 


by  line  500. 


TABLE  7.  PLOT/RESTART  TAPE  DATA  RETRIEVAL  PROGRAM 


00050 

00055 

00100 

00150 

00200 

00250 

00300 

00350 

00400 

00450 

00500 

00550 

00600 

00650 

00700 

00750 

00800 

00850 

00900 

00950 

01000 

01005 

01010 

01050 

01100 

01150 

01200 

01250 

01300 


' , LEN , NREC , REC ) 


/♦USERID  MVC01 
/♦TAPE  EDIT2 
//  EXEC  FGCG 
/ /SYS  IN  DD  ♦ 

LOGICAL* 1 CNE 
INTEGER*2  LEN 
LOGICAL* 1 REC ( 32000 ) 

NREC-  1 

CALL  PRECL('IN  ' , LEN , NR  EC , REC  ) 

DO  1 1-1,9999999 

IF(CNE(REC, 'PLOT' ,4)  )GO  TO  1 

KLEN-LEN/4 

CALL  OUT (REC , KLEN ) 

1 CALL  FBUF (499,488) 

99  STOP 
88  STOP 
END 

SUBROUTINE  OUT ( R EC , KL EN ) 

REAL  REC ( KLEN ) 

WR ITE( L 7, 10)REC(39) , REC (52) , REC (620), REC (621), REC (622), 
$REC (22  9) ,REC (1 169) ,REC ( 12  19) 

WRITE (6, 10)REC (39) , REC (52) ,REC (620) , REC (62 1 ) ,REC ( 622) , 
$REC (2  29) ,REC ( 1 169) ,REC ( 12  19) 

10  F0RMAT(F8. 3,  1 P7E  10. 3) 

RETURN 

END 

//DATA.FT17F001  DD  UN IT-B AT , FI L E S - $DEM 6 
//DATA-IN  DD  UN  IT- 2 4 00 , VO L- S E R- ED  I T 2 , DS N-T PO UT 
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When  a plot  record  is  found  program  execution  goes  to 
line  550  which  divides  the  record  length  LF.N  by  four, 
converting  the  record  length  from  computer  words  to  alpha- 
numeric words.  The  subroutine  OUT  Is  then  called  (line  600) 
and  execution  skips  to  line  850.  Line  900  changes  the  array 
REC  from  a logical  array  to  a real  number  array.  Lines  950 
and  1000  write  the  desired  parameters  (discussed  below)  to 
Fortran  Unit  17,  defined  In  line  1200  as  a BAT  file  called, 
in  this  case  $DEM6.  Lines  1005  and  1010  write  the  same 
parameters  to  Fortran  Unit  6,  a line  printer.  Following 
printing  of  the  desired  parameters,  program  execution  Is 
returned  to  the  main  program  at  line  650.  The  program 
continues  to  process  records  until  an  end-of-flle  Is 
encountered,  when  execution  Is  stopped. 

To  determine  the  word  number  of  the  desired  parameter 
to  be  printed,  the  following  steps  are  taken; 

For  system  parameters,  the  word  number  Is  simply  21 
plus  the  number  of  the  systen  parameter.  For  example,  the 
transient  time  Is  parameter  number  18  (recall  that  parameter 
numbers  are  taken  from  the  RELAP4  Users  Manual  on  the  pages 
noted  earlier).  Thus  the  word  number  for  transient  time  is 
21  + 18  ■ 39  and  REC(39)  contains  this  parameter. 

For  volume  parameters,  the  word  number  is  calculated  as 
f ol lows  , 


WN  - 41  (VPN  - 1 )NV0L  + VOLN 
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where  the  new  symbols  are, 

WN  * word  number  of  desired  parameter 
VPN  ■ volume  parameter  number 
VOLN  - volume  number  of  desired  parameter. 

For  example.  If  the  desired  parameter  Is  mixture  level 
in  volume  4,  the  word  number  Is  calculated  as  follows;  The 
VPN  for  mixture  level  la  9.  NVOL  - 23,  and  VOLN  - 4.  Thus 
WN  - 229,  and  RFC (229)  Is  the  mixture  level  In  volume  4. 

For  junction  parameters  the  word  number  Is  calculated 

f r om , 


WN  - 41  + 2 4 (N VOL  ) + (JPN  - 1)NJUN  + JUNN 

where  the  new  symbols  are, 

JPN  - Junction  parameter  number 
JUNN  • Junction  number  of  desired  parameter. 

For  example.  If  the  desired  parameter  Is  the  flow  rate 
in  Junction  27,  the  word  number  is  calculated  as  follows; 
The  JPN  for  flow  rate  la  1,  NVOL  - 23,  NJUN  - 33,  and  JUNN  - 
27.  Thus  WN  ■ 620  and  R EC ( 6 2 0 ) la  the  flow  rate  at  Junction 
27. 

For  heat  slab  parameters  the  word  number  la  calculated 

f rom  , 


WN  - 41  + 2 4 (NVOL ) +16(NJUN)  + (SLBPN  - 1)NSLR  ♦ SLRN 


where 


SLBPN  “ heat  slab  parameter  number 
SLBN  ■ heat  slab  number* 

For  example,  if  the  desired  parameter  Is  surface 
temperature,  heat  slab  S23,  the  word  number  is  calculated  as 
follows;  The  SLBPN  for  surface  temperature  Is  2,  NVOL  ■ 23, 
NJUN  - 33,  NSLB  - 25,  and  SLBN  - 23.  Thus  WN  - 1169,  and 
R EC (l  169)  is  the  surface  temperature  of  heat  slab  23. 


The 

central 

processing  unit  (CPU) 

time 

requl red 

for 

execut ion 

of  the 

program  depends  on  the 

amoun  t 

of  data 

to  be 

processed 

. As  a 

benchmark,  processing 

of  4 00 

records 

takes 

about  fifty  CPU  seconds. 


